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Abstract 

Random Matrix Theory (RMT) is the study of matrices with random variables 
determining the entries, and various additional symmetry conditions imposed on the 
matrices. A comparitively young theory, it has its roots in Hungarian physicist and 
mathematician Eugene Paul Wigner’s work in the early 1950’s [WigSla, WigSlb]. In 
the following years and continuing unabated to the present day, it has permeated nearly 
every area of modern physics and even number theory [Meh04, GMGW98, ABFll], 
In recent decades attempts have been made to further rehne what has become 
the canonical random matrix theory with its associated symmetry conditions, by 
considering symmetries which allow the matrix representations of quantum potentials 
to impose fc-body forces on the particles in a system containing m particles [k < m) 
[BRW01a,BRW01b,Sre02,Kotl4]. This is a generalisation of canonical RMT and is the 
topic of this thesis. It will be refered to as the unification of random matrix theories, 
because every randomised /c-body potential gives rise to a new ensemble. For k = m 
these ensembles are exactly those known and studied already under the rubrik of 
canonical RMT, but for fc < m the resulting ensembles are different and little is known 
about them. The phase space of random matrix ensembles generated by randomised 
/c-body potentials represents an opportunity to unify all random matrix theories into 
one, single mathematical theory of random matrices. The process involves embedding 
the ^-body potential into the m-particle state space creating what has become known 
as the embedded ensembles. The embedded ensembles, hrst introduced by Mon 
and French [MF75] in 1975, gave physicists a framework for studying many-body 
interactions using random matrix theory, but also presented practitioners with the 
challenge of an incomplete theory, as it showed that the (canonical) form of random 
matrix theory under consideration at the time was a single instance of a much larger 
phase space of random matrix theories. This meant that only a small subsection of 
the possible statistics were being calculated. Canonical RMT was already providing 
scientists, engineers and mathematicians insights into a variety of different patterns 
and physical phenomena in number theory, engineering, physics and computer science. 
However, these connections only became visible when they were uncovered, so the 
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potential that unification of random matrix theories offered presented both a big 
opportunity and new mathematical problems [BRW01a,BRW01b,Sre02]. Alternative 
mathematical methods would be needed to manage the new complexities manifested 
by studying the unihcation of random matrix theories. 

One of the greatest contributions to this area of study occured (in the author’s 
opinion) at the turn of the last decade with the publication of the breakthrough 
paper by Benet, Rupp and Weidenmiiller (BRW) which showed how a process of 
eigenvector expansions could be used to calculate certain statistical properties of 
fc-body potentials [BRWOla]. These methods were a great advance but were also 
difficult to implement, and it remains unclear if they can practically be used to 
calculate moments of the level density of embedded ensembles which are higher than 
the fourth moment. 

This thesis presents a completely different approach to the method proposed by 
(BRW). A new method involving particle diagrams will be introduced and developed 
into a rigourous framework for carrying out embedded random matrix calculations. 
Using particle diagrams and the attendent methodology including loop counting it 
becomes possible to calculate the fourth, sixth and eighth moments of embedded 
ensembles in a straightforward way. The method, which will be called the method 
of particle diagrams, proves useful hrstly by providing a means of classifying the 
components of moments into particle paths, or loops, and secondly by giving a 
simple algorithm for calculating the magnitude of combinatorial expressions prior to 
calculating them explicitly. By conhning calculations to the limit case m <C / —)■ oo 
this in many cases provides a sufficient excuse not to calculate certain terms at 
all, since it can be foretold using the method of particle diagrams that they will 
not survive in this asymptotic regime. Applying the method of particle diagrams 
washes out a great deal of the complexity intrinsic to the problem, with sufficient 
mathematical structure remaining to yield limiting statistics for the unihed phase 
space of random matrix theories. 

Finally, since the unihed form of random matrix theory is essentially the set of 
all randomised fc-body potentials, it should be no surprise that the early statistics 
calculated for the unihed random matrix theories in some instances resemble the 
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statistics currently being discovered for quantum spin hypergraphs and other ran¬ 
domised potentials on graphs [HMH05,ES14, KLW14]. This is just the beginning for 
studies into the held of unihed random matrix theories, or embedded ensembles, and 
the applicability of the method of particle diagrams to a wide range of questions as 
well as to the more exotic symmetry classes such as the symplectic ensembles, is still 
an area of open-ended research. 
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Chapter 1 


Introduction 


1.1 Physics and Unification 

Physics is the pursuit of natural laws, or models, which characterise and thereby predict 
the behaviour of events within our universe. This includes all measurable properties 
of the natural world, from how sub-atomic particles interact to the emergence of 
complicated molecules, proteins, plants and animals. Scientihc knowledge can be seen 
as a vast collection of these models and measurable properties, collected over centuries 
from the efforts, insights and gnesswork of all past scientists, both theoretical and 
experimental. Each scientihc discovery has its own unique history, interwoven into the 
fabric of previous discoveries. Sometimes a model discovered by a previous generation 
of scientists will be snperseded by a new theory. Many physicists acknowledge that 
the cycle may never end [FLS05,Boh60, Wei92,Bar91]. That it is possible that the 
pnrsnit of knowledge is a never ending process of destrnction and rediscovery, and 
that our understanding of natnre can always be improved, albeit never completed. 
The examples which epitomise this process of continual reformation is the replacement 
of Newtonian laws of motion by Relativistic laws, and separately the discovery that 
Newtonian Mechanics is in fact a limiting case of the mnch more general theory 
of Qnantum Mechanics [FLS05]. It also occurs that sometimes two models which 
previously were thought to be distinct are shown to be the resnlt of a larger more 
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general law. This is called unification. Unification is one of the most important 
aspects of scientihc labonr, since it reveals the connections between phenomena in 
the natural world which were previously thought of as being different, but which are 
in fact a result of the same underlying patterns. A seminal example of unihcation is 
the family of mathematical expressions called Maxwell’s equations^, which explain 
the relationship between electricity and magnetism, unifying two effects which were 
previously thought of as being distinct. 

Since the process of discovery often leads to different approaches to similar 
problems [Fey67] it is not always obvious that the underlying patterns governing 
different phenomena are the same. Not least of the problems is that seemingly 
disparate scientihc questions are often studied within the conhnes of their own specihc 
jargon. To break down the barriers between problems and to reveal their in h erent 
similarities it is essential to look uncompromisingly at the underlying mathematics 
without being swayed by the layer of human language which has been added “on top”. 
Another essential feature of unihcation theories is that they include mathematical 
parameters which connect things which were previously thought of as being unrelated. 
In this sense unihcation requires an expansion of the mathematical expressions to 
“make room” for effects which were not given a voice in previous theories. This thesis 
will look at a particular case of unhnished unihcation in physics and mathematics; 
the unihcation of random matrix theories. It will do this hrstly by observing that 
the common, or canonical, form of random matrix theory is just a single point in 
a whole landscape of possible random matrix theories. This landscape of random 
matrix theories is known as the embedded ensembles (see Chapter 2) which are 
themselves a subclass of quantum fc-body potentials. In Chapter 3 it will be shown 
heuristically how Wigner’s Semi-Circle law for canonical random matrices can be 
calculated using diagrammatic methods. In the same chapter an illustration of the 
method of supersymmetry will show the sometimes stark difference in mathematical 
approaches to the same problem, in this case Wigner’s Semi-Circle Law. Finally in 

fii) V • E = ^ (ii) V • B = 0 (iii) V x E = (iv) V x B = //q (J + eo^)- Together 
equations (iii) and (iv) describe the electric field as a function of the magnetic field and vice versa, 
showing that E and B are two aspects of a single unified physical process. 
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Figure 1.1; The canonical form of random matrix theory is just a single point in a 
whole landscape of possible random matrix theories (unified RMT). Unified RMT in 
turn represents the set of randomised quantum many-body potentials, also known as 
the embedded ensembles. 


Chapters 4 and 5 it will be shown that a method of particle diagrams, involving 
graphs of the relationships between quantum states, is sufficient for calculating the 
moments of the level density of the embedded (unified) ensembles. The majority of 
the results presented here have been published in the rapid communication [SM14] 
and the article [SM15]. 











Chapter 2 


Random Matrix Theories 

2.1 Canonical RMT 

2.1.1 Statistical Landscapes 

In the early 1960’s theoretical physicist Freeman Dyson pnblished a trio of papers 
with the heading Statistical Theory of the Energy Levels of Complex Systems [Dys62a, 
Dys62b, Dys62c] . The results (which were also inspired by earlier works of Eugene 
Wigner and Madan Lai Mehta [WigSla, WigSlb, Meh04]) and the related held of 
mathematics which grew around them, became known as random matrix theory 
(RMT). The basic tenet of random matrix theory as introduced by Dyson was the 
investigation of the statistics of “a// physical systems with equal probability". To 
paraphrase events, a challenge which researchers were facing at the time was to 
characterise the behaviour of large complex nuclei. Wigner, Dyson and Mehta 
envisioned approaching the problem by assuming nothing, or at least as close to 
nothing as possible. From this position they began to investigate the quantities that 
can be known about the behaviour of a particular physical system. Their starting 
point was the hamiltonian of the physical system with the single constraint 

H = H^ ( 2 . 1 . 1 ) 
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for physicality. Assuming nothing else, one would not be able to say anything about 
a particular hamiltonian H. However, by allowing each H to occur with some 
probability p{H) it becomes possible to study the whole space of probabilistically 
weighted hamiltonians. From this point onwards the investigation necessarily becomes 
statistical, and these ensembles of hamiltonians are stndied by asking qnestions of 
the form “What is the probability that... (etc. )^\ The proposition which then needs to 
be verihed with experiment is that a natnral system is some nnknown hamiltonian 
obeying (2.1.1) and therefore the typical behavionr of the system may obey the 
behaviour specihed by the ensemble average of this. In other words, althongh the 
information extracted from these ensembles of hamiltonians takes the form “ What is 
the probability that...^\ it is possible to use the probabilistic results to say something 
about the properties of the “typical” (read average) of a hamiltonian, because 

f{X) = J f{X) ■ (the probability that X) dX (2-1-2) 

where the horizontal line above f{x) indicates the average taken over the ensemble of 
X’s. Investigating random hamiltonians in this way is what gave fruit to the field of 
random matrix theory - a family of sophisticated mathematical tools and equations 
for extracting as mnch information as possible abont an ensemble of hamiltonians 
starting from only a handful of very basic assumptions [Meh04]. The additional 
assumptions would usually take the form 

1. H = 

2. H = S-^HS (2.1.3) 

3. Gaussian probability determining the matrix elements of H 

where the matrix S is some predehned matrix determined by the particular physical 
setup at hand. It imposes some symmetry on H additional to the hermitian symmetry 
it already satishes by virtue of (2.1.1). The three classical examples for S as introduced 
by Dyson will be given in the next section (one of which simply takes S = 1). These 
essentially divide the physicial landscape into four statistical landscapes; three groups 
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defining a “threefold way” and a fourth being the complement of these (everything else 
not included in the three groups determing the threefold way). Although the initial 
development of a theory of random matrices focused on the first three symmetry 
groups defined by Dyson, it turns out that this classification defines just a single 
point in a much larger landscape of random matrix theories and a more nuanced 
approach to symmetry is required to unify them all. 

2.1.2 Symmetry Classes 

Dyson’s initial development of random matrix theory was placed within the context 
of a “threefold way”. These are three symmetry classes which determine a set of 
the additional symmetries (see (2.1.3)) obeyed by physical hamiltonians. Recreating 
some of Dyson’s steps, this section will detail the formulation of the three symmetry 
classes determining the additional constraints on the matrix elements of H when the 
hamiltonian refers to a quantum system from one of three specific groups. Later it 
will be shown why this analysis is incomplete, constituting just a single point in a 
phase space of additional symmetries. The three groups are distilled from the set of 
all possible hamiltonians by asking some fundamental questions. Eor example it is 
reasonable to ask the question “What form does the potential take if the physical 
system contains only fermions, or only bosons?” Similarly, it is natural to wonder 
“What form does the potential take if the system is time-reversal invariant compared 
to the case when it is not?” These are questions of symmetry which are integral to the 
original purpose of RMT, which is to assess the statistical properties of ensembles of 
hamiltonians. Although random matrices do not have to be hamiltonians per se - they 
can also be scattering matrices [ABEll], or any other matrix for that matter - the 
classical categorisation of RMT into three symmetry classes as given by Dyson takes 
its starting point as the hamiltonian of a physically permitted quantum system. In 
other words, the starting point is the hamiltonian satisfying the Schrbdinger equation 
of the system 

dih 

Hip = ih— = Eip. 

ot 


(2.1.4) 
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Even with a modest amount of complexity Schrodinger’s equation becomes too 
difficult to solve analytically, which is why RMT was developed in the first place. 
In RMT it is assumed that the Hamiltonian H is highly complex and only some 
rudimentary symmetry properties of the physical system are known. After placing 
these restrictions on the hamiltonian all remaining free parameters are randomised. 
Thereafter an attempt is made to calculate whatever statistical properties can be 
yielded as a result of the randomisation and symmetrisation process. Using this 
probabalistic view the problem becomes one of studying ensembles of matrices, each 
grouped into families depending on the set of minimum conditions imposed upon them 
and each weighted by some probability of occurring within the ensemble. Instead 
of studying the eigenvalues of a single instance of H one calculates distributions of 
eigenvalues defined by how likely they are to occur across the entire ensemble. 

Particle Wave Functions 

The starting point for calculating the three symmetry classes proposed by Dyson is a 
remarkable law of nature which through repeated experimentation has established 
that any quantum state {a.k.a wave function or p.d.f) describing a physical experiment 
containing elementary particles is either even-symmetric or odd-symmetric under 
particle exchange. Particles with even wave functions (bosons) and particles with 
odd wave functions (fermions) can be represented mathematically with the respective 
expressions 

^7(rpi,... ,rp„) = ij{ri, ... ,r^) (bosons) 

^/>(rpi,...,rp^) = (-l)^^/>(ri,...,r,„) (fermions) (2.1.5) 

where {Pi,..., Pm} is a permutation of (1,..., m} and P is the number of pairwise 
permutations between the single particle states required to bring (Pi,..., Pm} back 
to the initial configuration. More concisely, for the wave function of both bosons and 
fermions the following symmetry condition holds 


fj{rpm...,rpJ = ^^fj{ri,...,rm) 


(2.1.6) 
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^Bosons = 1) 

Fermions = — 1) ^ 

photons 

protons 

pions 

nentrons 

mesons 

electrons 

gluons 

muons 

^He 

nentrinos 


quarks 

V 

■^He J 


Table 2.1: Table of common bosons and fermions, the two classes of quantum particles 
in nature. Bosons can always be represented by an even probability density function 
while fermions, obeying the Pauli Exclusion Principle, can always be represented by 
an odd probability density function. 


where as before = 1 bosons and — 1 for fermions. A composite particle 
composed of any nnmber of bosons and an even(odd) nnmber of fermions will 
have an even(odd) state representation. That is, the composite particle itself be 
a boson(fermion). Hence it is possible to divide np the vast physical landscape of 
possibilities into just four groups; 

1. Time Reversal Invariant Bosonic systems. 

2. Non Time Reversal Invariant Bosonic and Fermionic systems. 

3. Time Reversal Invariant Fermionic systems. 

4. None of these. 


The first trio of systems define the three symmetry classes of canonical RMT. One 
may protest that all possible realities are taken into account by the hrst three groups. 
However, each group is defined in terms of the minimum set of restrictions on each 
member so systems reqniring additional restrictions mnst go into the fonrth gronp. 
What are the attendant symmetries that mnst be satished by a potential H for a 
system in each gronp? 
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Invariance under Time Reversal 

To form a meaningful mathematical picture of how time reversal invariance relates 
to a physical measurement it is first necessary to look at unitary operators, or more 
precisely in the present context, unitary matrices. A good introduction to the following 
ideas can also be found in [Meh04, HaalO]. A unitary matrix U, is a matrix which 
conserves the overlap between states 

WI170) = im- (2.1.7) 


Rearrangment gives 


UW = 1 . 


( 2 . 1 . 8 ) 


Expressed in terms of a summation over the individual elements of the matrix this is 
equivalent to 

T.Ua,u;, = s„. (2.1.9) 

k 

In other words the rows (columns) of a unitary matrix are orthonormal. Replacing 
(f) with "0 in (2.1.7) gives ||([/V’|tlV’)ll = 11 ('01V’) 11 so that Unitary transforms also 
conserve the length of a vector. From the above it follows that U G rigidly 

rotates a vector in C” while conserving its length. Notably, a state transformed as 
0 —>• U'ljj is indistinguishable from the original state when the measurement operator 
is also transformed as 

A UAt/0 (2.1.10) 

Additionally assuming a system of particles whose states are of the form 'ip = 'll)^ 
it follows immediately that the overlap between two arbitrary quantum states is 




( 2 . 1 . 11 ) 


If time is flowing in the forward direction the coordinates of a particle will be given 
by its position, momentum and t-variable {x,p,t} whereas if time runs backwards 
the particle’s momentum and time variable will reverse at every position so that the 
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coordinates become {*, —p, —t} and the overlap between two states will be 

= im* ( 2 . 1 . 12 ) 

where \'ip^) is the wave function l"^) under time reversal. Representing the time 
reversal of a state as an operator gives the tautological expression 

Ti! = (2.1.13) 

where by (2.1.12) T must obey 

(r^|T0) = {^\ct>Y. (2.1.14) 

That is, the operator T is anti-unitary. Additionally, under time reversal the Hamil¬ 

tonian of the system becomes 

rR ^ tHT-\ (2.1.15) 

Comparing (2.1.14) with (2.1.7) it can be seen that the only difference is complex 
conjugation, so T can be written as 


T = KC (2.1.16) 

where K is unitary and the operator C is defined by its action of taking the complex 
conjugate of the state on which it operates. Since double application of time reversal 
to a state should leave it in a state indistinguishable from itself i. e. differing only in 
phase, T must without exception obey the relation 

= a (2.1.17) 


with |a| = 1, so 


KK* = a. 


(2.1.18) 
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Additionally, by (2.1.10) T transforms as 

T ^ UTU-^ (2.1.19) 

under unitary transform of the wave-function by U, and the matrix K dehned 
previously transforms as 

K UKCU-^ = UKC {U^y = UKU^. (2.1.20) 

A system is called invariant under time reversal if 

H^ = H (2.1.21) 

where H is the hamiltonian. A matrix satisfying (2.1.21) is called self-dual. A time 
reversal invariant system therefore satisfies the following sequence of equalities 

H = THT-^ = KCH{KC)-^ 

= KCHCR-^ = KH^K-^ ( 2 . 1 . 22 ) 

where the fact that H is hermitian has been used, so that H* = . Recalling that 

K is unitary, implying K*K^ = 1, and motivated by a stroke of insight to multiply 
(2.1.18) by it is seen that 

K = aK^ = a^K (2.1.23) 

which means that the value of a is restricted to satisfying 


a = ± 1 . 


(2.1.24) 


Twofold Way 

The insight that a can take two possible values reveals mathematically that there 
are two distinct types of systems, both satisfying time-reversal invariance. Equation 
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(2.1.24) implies two possible outcomes 

KK* = 1 or (2.1.25) 

KK* = -1. (2.1.26) 

In addition to being nnitary, if (2.1.25) holds K is symmetric and (2.1.26) implies 
K is anti-symmetric. Eollowing, bnt not proven here, is that in the symmetric case 
there is a nnitary transform which by (2.1.20) gives 


K = 1 


(2.1.27) 


which implies that the system is bosonic (even spin^) and in the anti-symmetric case 
there is a nnitary transform snch that K becomes 


K 


/ 0 +1 

-1 0 

O 

V 



0 +1 
-1 oj 


(2.1.28) 


which implies that the system is fermionic (odd spin^) [Meh04,Haal0]. That is, a basis 
can be chosen -0 —)■ U'll) for each of the two possible time reversal invariant systems 
snch that K has the above forms for each case respectively. After choosing snch 
a basis snbseqnent transforms on K can only be performed if the relevant relation 
(2.1.27) or (2.1.28) remains trne. Hence for bosonic systems fnrther transforms are 
restricted by (2.1.20) and (2.1.27) to being of the form 


K ^ OKO'^ 


(2.1.29) 


^This thesis will use bosonic and fermionic interchangeably with, respectively, even spin and 
odd spin throughout. 
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where O must be orthogonal (real unitary) and hence by (2.1.22) H is real hermitian 
symmetric. To summarise, the above developments gives justihcation for the following 
definition 

Definition 1. The ensemble of Hamiltonians of Bosonic time reversal invariant 
systems are hermitian symmetric matrices statistically invariant under 

H OHO'^ (2.1.30) 


where O is orthogonal. 

Here statistical invariance means that the probability of a given H occuring within 
the ensemble is the same as the probability of any other H found after applying the 
transform, which in this case is (2.1.30). The ensemble dehned by Dehnition 1 along 
with the condition that the p.d.f. of the elements of a member of the ensemble are 
statistically independent dehnes the Gaussian Orthogonal Ensemble (CUE) which is 
denoted by /9 = 1. 

Turning now to the fermionic case note that by (2.1.28) the number of rows 
(columns) of K is always even. Once a basis is chosen such that (2.1.28) holds, further 
transforms on K are restricted to satisfying 

K = ZKZ'^ (2.1.31) 

where Z is unitary. Matrices Z satisfying (2.1.31) form what in the literature is called 
the symplectic group. Notice that by substituting (2.1.28) into (2.1.22) it can be seen 
that the hamiltonian of a time reversal invariant Eermionic system satisfies 

(THT-\, = = H., = Hi, 


(2.1.32) 



14 


CHAPTER 2. RANDOM MATRIX THEORIES 


which can be written as 


0 

Oil 

bn 

^12 

bi 2 ■ ■ 

® im 

blm 

ei 

—h* 


—h* 

^12 

®12 

■ —h* 



62 

0 

0-21 

621 ■ ■ 

®2m 

^2m 



62 

—h* 

^21 

a* 2 i ■■ 

—h* 

^ 2 m 

®2m 


H = 


* 


V 


0 

6m j 


(2.1.33) 


where the * in the lower diagonal denotes that the matrix is hermitian so that the 
lower diagonal is defined by the upper diagonal elements. The diagonal of H is made 
up of 2 X 2 blocks of the form 


e 0 
0 e 


(2.1.34) 


where by hermitian symmetry e must be real. The off diagonal components consist 
of 2 X 2 blocks of the form 




(2.1.35) 


where a, 6 G C. A useful shortcut to the same conclusion is to calculate the matrix 
elements of H relative to the basis 


T'iIj 2 , ..., T'lp2n- (2.1.36) 

Since T = KC and for this (symplectic) case = — 1 it follows straightforwardly 
that 


(^|iP|T^) = - WITTI^) 


(2.1.37) 

(2.1.38) 
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which is the equivalent of (2.1.35). Taking summary once again, the above discussion 
gives motivation for the following definition. 

Definition 2. The ensemble of Hamiltonians of Fermionic time reversal invariant 
systems are hermitian self-dual (H^ = H) matrices statistically invariant under 

H ZHZ-^ = ZH{KZ^K-^) = ZHZ^ (2.1.39) 

where Z is symplectic. 

Matrices which additionally have the property that entries not related by symmetry 
are statistically independent form the Gaussian Symplectic Ensemble (GSE), which 
is denoted by /? = 4. 

Non Time-Reversal Invariance 

The Hamiltonian of non time-reversal invariant Bosonic and Eermionic systems is 
unrestricted other than by statistical invariance under a unitary transformation. 
Although this has been a guiding restriction for the cases /? = 1 and /9 = 4, the 
additional demands of time reversal invariance confined these to having the symmetries 
defined above. Relaxing these restrictions, which implies S' = 1 in (2.1.3), gives the 
following definition of non time-reversal invariant systems. 

Definition 3. The ensemble of Hamiltonians of non time-reversal invariant Bosonic 
or Fermionic systems are hermitian matrices statistically invariant under 

H UHU-^ (2.1.40) 


where U is unitary. 

If in addition the entries of the matrices not related by symmetry are statistically 
independent Definition 3 defines the Gaussian Unitary Ensemble (GUE), which is 
denoted hj (d = 2. Definitions 1, 2 and 3 form Dyson’s threefold way. They give 
the restrictions on the matrix representation of the potential H depending on which 
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class of quantum system it describes. By classifying the physical landscape into these 
classes the matrix H is found to obey additional symmetry (on top of if = ii^), even 
if everything else remains, for the present, unknown. The symmetrising matrices O, 
Z and U dehne the matrix S of (2.1.3). It should be emphasized that the additional 
symmetries satished by H in each of these three groups are the minimal constraints 
satished by any H belonging to the set. Hence the set of non time-reversal invariant 
Hamiltonians dehned by the GUE do not exclude the time-reversal invariant systems 
dehned by the GOE and GSE ensembles. Those matrices are present in the ensemble 
of GUE matrices as well, but occur with a lower probability. The set of non time- 
reversal invariant hamiltonians is then more strictly the set of “not necessarily but 
possibly time-reversal invariant” hamiltonians. 

“None of These” — Embedded RMT 

While the canonical approach to random matrix theory places a single random variable 
in each cell of the matrix, it will be shown in subsequent sections that there are in fact 
additional ways to dehne the potential H, these being dehned in terms of the order 
of the potential, k, which will be dehned later under the framework of many-body 
potentials. In this way the canonical form of RMT will be extended and a new model 
will be proposed which has the canonical form as a special case. Under the new model 
what is normally referred to as RMT will be shown to be just one of the possible 
points in a phase space of random matrix theories, each with its own set of statistical 
properties. This unihed phase space is sometimes referred to as Embedded Random 
Matrix Theory. 

Instead of looking at classical, or canonical RMT hamiltonians, the purpose of 
this thesis will be to investigate the statistics of Embedded RMT hamiltonians. These 
dehne a unihed form of RMT which is still in its infancy; each embedded hamiltonian 
represents a particular instance of a random matrix theory, and a unihed theory 
of random matrices must classify the statistical properties of each one. This thesis 
will show how the symmetry restrictions imposed in any particular instance of a 
random matrix theory can be represented in terms of particle diagrams and a new 
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mathematical methodology will be introduced to make calculations of the statistics 
of quantum many-body systems from a range of random matrix theories in the phase 
space. 

2.1.3 Many-Body Potentials 

It will be seen that the hermitian matrices of canonical RMT with GOE, GSE or GUE 
symmetry imposed on the matrix elements of the hamiltonian H represent only one 
possible flavour of a random matrix theory. To address this problem and in an attempt 
to unify the held of random matrix theories, K. K. Mon and J. B. French introduced 
the embedded RMT ensembles [MF75], which are hamiltonians written in second- 
quantised form and determined by a trio of parameters k, m and I (dehned later). 
For each set of values {k,m,l} one attains a distinct random matrix theory. These 
second-quantised hamiltonians were already studied by physicists before the advent of 
RMT in the context of Many-Body Quantum Mechanics where they appeared under 
the rubric of Many-Body Potentials. In this section the intention is to describe in some 
detail what a quantum many-body potential is, which will involve an introduction 
to the basic notational norms of second-quantisation. A good induction into the 
following formalism can also be found in [NO88,DN05]. The preliminary aim will be 
to express an arbitrary many-body operator in terms of creation and annihilation 
operators. The hamiltonian operator H expressed in this way will form the extension 
to canonical RMT under consideration. It will then be seen that there is a particular 
set of values for the parameters {k,m,l} determining this hamiltonian where it returns 
to the form used in classical RMT which, as noted, is a specihc case of the more 
general class. 

In the forthcoming model the number of quantum particles in the system of 
interest is a variable denoted by the letter m. The tensor product of m single-particle 
states 

|q!i ... a™) = |ai) ® ® |ttm) (2.1.41) 

gives a natural way of describing an m-particle state. A simple list of the single¬ 
particle states is written in a single packet {ket) and dehned as the quantum state. 
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In position representation this m-particle state can be rewritten as 

= {ri,...,rm\ai...am) = (o|ai)(r2|a2) • • • (a^lr^) 

= 0ai(o)0a2(»^2) • • •0a^(rm)- (2.1.42) 

Given an m-particle state ipm the properly symmetrized and normalised bosonic {B) 
and fermionic (F) state is therefore given by 

P{B,FMri,...,rm) = (2.1.43) 

where the symbol .^ = 1 for states describing bosons and — 1 for states describing 
fermions. P is the set of permutations on {1,2,..., m} and the same letter is also 
used to signify the parity of the permutation. To simplify notation later, dehne the 
properly symmetrized (but not necessarily normalised) m-particle state by 

|ai. ..a^} = \^.P{B,F}\ai ■ • • a^) = (2.1.44) 

Vm- p 

The closure relation for this space is then given as 

'y ) \P{B,F} |cii . . . Clm)] [P{B,F} | Cli . . . Clm)] j" ^ ) | ttl . . . CUrn} {oil ■ ■ ■ Oim \ 

= 1. (2.1.45) 

Since m-particle states are regarded as describing a probability density function they 
should be normalized to unity. The odd symmetry of fermionic states prohibits 
them from containing any two particles with the same state, so the overlap of two 
m-particle Fermionic states when non-zero is 

(a) a^} = (-1)^- (2.1.46) 

If the two states do not contain the same number of particles there is certainly a zero 
overlap. Bosonic m-particle states on the other hand - having even symmetry - are 
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allowed to have single particle states equal to other single particle states. This means 
that when the overlap is taken between two bosonic states containing k unique states 
with Hi single particle states in state i,l < i < k the result is 

{a\ ... aj\ai... = ni\n 2 \... n^! (2.1.47) 

where the number of particles is m = One can fortunately condense this 

notation further, expressing (2.1.46) and (2.1.47) as the single expression 

{a[ ... a'^\ai... am} = (2.1.48) 

a 

Hence for functions representing a collection of bosons the orthonormal states are 
given by 

1^1 ... am') / ^ I• • • O^m} (2.1.49) 

VHa nj 

whereas for states representing a collection of fermions this simplifies to 

\a \... am) \^i • • • ^m}‘ (2.1.50) 


To conclude, the orthonormal closure relation for bosons and fermions is 


I ^ ^ I®1 • • • ^m} {oil ■ ■ ■ Olm\ ^ ) 

ml 


ria nc 


ml 


Icr 1... am) • • • ^m\ 1* (2.1.51) 


fc-Body Operators 

The above formulates a way to describe states which represent collections of bosons 
and fermions. How to mathematically describe the Hamiltonians which exert forces 
on these states? The hamiltonians will be expressed as A;—body operators, where k is 
a variable determined by the particular system under consideration. An operator is 
called a /c-body operator if its effect on an m-particle state is the sum of its effect 
on each of the single particle A:-tuples contained in the m-particle state. For 
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example a 1-body operator is an operator satisfying 

m 

H\oii... Oi-m) 'y ' ... OijYi') (2.1.52) 

i=l 

where Hi is the effect of H hitting only the i’th single-particle state. It can therefore 
be seen that 


(tti... ttffjI Tr|/3i ... /3j^) y ) (tti... I fPj|/3i ... /3m) y ) {o^klPk}{p^i I l/3i) 

i i k^i 


(2.1.53) 


and for non-orthogonal \ai... am) and |/3i... (3m) one can conclude that 


(ttl ... am\H\(3\... /3m) _ {oii\H\/3i) 


(tti . . . am\(3l • • • Pm) 


{oiilPi) 


(2.1.54) 


The dehnition for 2-body operators is the natural extention of this idea. The effect of 
H 2 on an m-particle state is the sum of its effect on each of the single-particle 
pairs contained within the composite m-particle state so that 

H\ai... am) y ) 33^1 j |cii... am) 2 ^ ' 33ij |cii... am/)■ (2.1.55) 

i<j ^ 

As with the 1-body case one can write the following relation 

(Ol . . . am\H\P\ . . . Pm) X y ( (o^l • • • 0!m\Hij\Pi . . . Pm) 

^ Hi 

= \'y2I[ {ak\Pk){aiaj\Hij\PiPj) (2.1.56) 


and diving across by non-orthogonal |ai... |/9i • • • /9m) gives the expression 


... cXyYi\H\l3i... /3j7^) 
(q;1 . . . 0(rn\Pl • • • Pm) 


1 V (^0(iCXj\Hij\PiPj) 
2^ {ai\Pi){aj\Pj) ■ 


(2.1.57) 
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Following in direct analogy to the 1- and 2-body potentials above, the dehnition for 
the general /c-body case is given by 


H I CTi . . . OiffiJ 'y ' Lfjj jj, I Q?! . . . Cy-rn) 


(2.1.58) 


where the result of hitting the entire m-particle state means the same thing 

as H hitting the fc-body state \ai^ ... on^)- The normalized matrix elements of the 
operator are therefore given by 


(q?! . . . (y.jYi\H\(3\ . . . firn) _ 1 

(Oi . . . (y.rn\l^l • • • l^m) 


(Q^U • • • 0^4 1/^4 • • • /^ik) 

(«4lAi)---(a4lAJ 


E 


4^*2...^4 


The above explains how the many-body potential H acts on many-body quantum 
states, but what is H actually, mathematically speaking? 

Second Quantisation 

To describe H mathematically it is necessary to reformulate its properties as described 
already (which is to say, how it acts) in terms of creation and annihilation operators. 
These form an efficient way of talking about m-particle states because by using 
creation and annihilation operators it subsequently becomes possible to talk about 
the existence or non-existence of the single-particle states comprising each m-particle 
state. This binary way of thinking matters. It opens up a variety of shortcuts, for 
example when taking overlaps between m-body states; the answer can be calculated 
simply by looking at whether certain single-particle states exist in the m-body states. 
This “superpower” will prove extremely useful later on. 

To begin, the m-particle creation operator a\ adds a particle in the state |A) to 
the m-particle state on which it acts 



(2.1.60) 


This is consistent even if A is already present in the m-particle state. By (2.1.49) it 
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follows that 


giving 


®aI«i • • • «m) = 'Jnx + l^Y[nJ\Xai... am) 


• • • O^m) \/^A T IjAcii . . . am) 


(2.1.61) 


(2.1.62) 


where nx is the number of single-particle states equal to |A) in the original m-body 
state \ai.. .am)- Dehning |0) as the vacuum state, any other state can then be 
expressed in the form 




“Llo>- 


(2.1.63) 


To express an arbitrary fc-body operator in terms of creation and annihilation operators 
in the standard form the commutation relations between them are needed. Eirst it 
can be noticed that since 


^A^ul"^! ■ ■ ■ |A/iAi... An^j” ^IpAAi... A^} ^oj^ci^lAi... (2.1.64) 


one 


has which yields 


[«A>«i]-? := 44-^44 = 0. 


(2.1.65) 


Taking the complex conjugate gives 


[®A) ^ 0- 


( 2 . 1 . 66 ) 


Equations (2.1.66) and (2.1.66) give the first two commutation relations. To hnd the 
third requires an investigation into the operator ax '■= (al)^- Taking the overlap 


{cii . . . |/3i . . . fin\ {Aoi . . . Clm|/3l . . . (^n\ ■ 


(2.1.67) 
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one sees this can be non-zero only if m -|- 1 = n. Representing the action of ax in 
terms of the identity over all possible quantum states i.e. including states containing 
different numbers of particles, and assuming that |A) is in |/9i... /9„}, gives 

CX) 2 

^a|/3i • • • Pn} ^ ^ j" ^ ^ • • • /3rx}'|Q;i . . . CAm}' 

CO 

^ ^ r ^ ^ . . . cXfyi\/3i .. . 

777 ' 

m=0 <71...Qm 

n 

i 

n 

= i:r"i5AS.IA...A...W (2.1.68) 


where /3j denotes that the i’th particle is removed from the set. Hence the effect of 
a\ on the state is to remove a particle in the state A. The operator ax is therefore 
referred to as an annihilation operator, its normalized formulation being 

aM ... A.) = -jj= y I ,/nT I 

i VV i / 

1 ” 

= ^ y r-'^Ae.l/3i... ft ... /3„>. (2.1.69) 

Y ^A I 

As an aside note that for bosons (.^ = 1) this becomes just 

aA|n/3i ... ua ... np^) = ... (ua - 1)... np^) (2.1.70) 

where ... np^) denotes the state with np^ particles in the state jSi. It should be 
noted that the hnal commutation relation [oa, is still unknown. It can be found 
hrstly by noticing that 

axa\\ai... a„} = ax\irai... a„} 

n 

• • • C^n} ^ ^ ^Aaj . . . CXi . . . CXn} 

i=l 


(2.1.71) 



24 


CHAPTER 2. RANDOM MATRIX THEORIES 


and likewise that 

n 

* * * ^n} ^ ^ | q ^1 • • • (y.i . . . Otn\ 

i=l 
n 

= ... aj. (2.1.72) 

i=\ 

Combining these two equations gives 

^A/j-|*^l • • • ^n} H“ Clj^Cl_x|Ofi . . . Cy-n} 

= [5a/. + ^ala\] |ai • • • On} (2.1.73) 

which is the final commntation relation songht. Namely 

[ax,alU = 6,,. (2.1.74) 

Given (2.1.65), (2.1.66) and (2.1.74) an expression can now be found for an arbitrary 
fc-body operator in terms of creation and annihilation operators, initially by working 
in a basis where the operator is in its diagonal form, and hnally by expressing the 
result in any basis. To begin, following (2.1.58) and assuming that Hk is diagonal 

(c^ii • • • Cljj, I Tffc I Q!jj ...CTjj.) (2.1.75) 

SO that for arbitrary m-body states \ai^ ... a*^) and l/dj^... A^) the following sequence 
of equalities are attained 

m 

(ttjj . . . . ./3j^) — (ttjj . .. Q!j^ I — ^(0 .../3j^) 

= n ('^il/^i)('^n • • •'^*fel-^*i...*felAi • • • Afe) 

~ ^\ n I A) (Ai ■ ■ ■ Afe l-^fcl Ai • • • Afc) (®ii ■ ■ • *^41 Ai • • ■ Afc) 

— "^(cTil • • • Q!im I Ai • • • Am) (Ai • • • Afe l-^fcl Ai • • • Afc) 


(2.1.76) 
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the sum being over all ^-tuples in the state (3 := \(3i^ ... 

Tuple Counting 

Expressing (2.1.76) in terms of the matrix elements of the operator gives 

(®il • • • ■ ■ ■ Am) ~ ^ 1 ■ • I Al ' ' ' An ) (2.1.77) 

ku—tuples 

the sum now running over unique fc-tuples (any permutation of a given /s-tuple is 
considered non-unique) and the tuple-counting operator giving the total number 
of fc-tuples denoted by the labels fc = {A ... A} in the state upon which it acts. For a 
1-body operator the number of particles in the state a is in fact given by the number 
operator Ua = so that 

Ti = a^aa (2.1.78) 

and similarly the number of pairs consisting of the single-particle states |a) and |/3) 
is n^np if a A and nain^ — 1) otherwise. These can be combined into the single 
condition 

T 2 = na(na - Sap) = aialapUa (2.1.79) 

Generalizing immediately for a fc-body operator the number of fc-tuples are given by 

the operation 

Tk 'y ) (?7.j2 *^*1*3 ^*2*3) • • • (^*fe *^*2*fc ■■■ ^ik-li-k)' (2.1.80) 

Using the identity 

^pSpq (^qSpq (2.1.81) 

equation (2.1.80) reduces neatly to 

Tk — 


( 2 . 1 . 82 ) 
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so that the second-quantized formulation of a diagonal fc-body operator is 

H = ^ Y. • • • «n- (2.1.83) 

* Z1...4 

To express this in terms of an arbitrary basis observe that in terms of some basis |A) 
one can write 

(2.1.84) 

a 

Taking the complex conjugate gives 


= 5I(«Aja)( 


(2.1.85) 


and substituting these expressions into (2.1.83) leads to the following 


H ^ Y ^<...4 E0'iIh)«5i • • • Y(jk\i'k)a], Yi^Mai,... 






3k 


'^k 


~ 'Yj 'Yj • •-ifclH• • ■ r)(h• • • 


6 ■ 6 7 
1 fc ll 




• • • ^jk^ik • • • 


( 2 . 1 . 86 ) 


n---ok 

h--ik 


which is the second-quantized form of a fc-body operator in an arbitrary basis. 
This formula expresses the hamiltonian matrix H of (2.1.1) as a second-quantised 
operator i.e. an operator defined by a sequence of creation and annihilation operators. 
Implicit in the model are three important parameters. Firstly there is m, the variable 
determining the number of single particles in the system. In other words for any state 
\fi) one has 

\fi) = \ai... am) (2.1.87) 

for some set of single-particle states {«!,..., am} as seen in (2.1.50). Secondly there 
is k < m, the order of the interaction or in other words the “number of bodies” 
involved in each interaction under the force of the potential. Einally, there is the 
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implicit parameter I which determines the number of energy levels available to each of 
the single-particle states in the compound m-body state. Hence I is the size of the set 
from which the single-particle states ai,..., can take their values. As mentioned 
in section 2.1.3 the trio of values together with the symmetry conditions 

imposed on H form a single instance of a random matrix theory. As an unrestricted 
phase space the hamiltonian of (2.1.86) represents the unification of these random 
matrix theories. Next it will be shown that canonical RMT coincides with the case 
k = m with I — )■ cx). 

2.1.4 Canonical RMT as a Single Point 

To paraphrase section 2.1.1 Canonical RMT is the study of random matrices with 
some predefined symmetry conditions and a single random variable determining 
the quantity in each cell of the matrix. The unified form of random matrix theory 
determined by the hamiltonian of (2.1.86) however, allows for the possibility that 
more than one p.d.f. determines any given cell of the matrix. For the special case 
where k = m this becomes 

^ ~ TiT ^ (2.1.88) 

SO that for any two m-body states 

1/^) 1*^1: * * * : ^m) (2.1.89) 

= (2.1.90) 

the matrix elements of H become 

Hfii, = y | ^) 

il .Jm 


(2.1.91) 
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Hence the cells of the matrix H for the special case k = m contain only one element 
- the hnal line of (2.1.91). There is no summation as there would be for the case 
k < m. This single element, being dehned in terms of a probability density function 
and symmetrised by some condition on H will of course give back the corresponding 
canonical RMT ensemble symmetrised under the same condition. This tells us that 
k = m IS the point in the phase space {fc, m, /} of the unihed theory which coincides 
with canonical random matrix theory. Moreover, since the random matrices of the 
canonical theory are usually assumed to be inhnite one takes / —)■ cx), giving canonical 
RMT as the theory coinciding with the point {m, m, cx)} in the unihed phase space. 




Chapter 3 


Wigner’s Semi-Circle Law 


3.1 Wigner’s Law for k = m 

Wigner’s Semi-Circle Law is one of the more iconic and widely known resnlts to come 
out of the held of random matrix theory. This is the rule, proven mathematically, 
which states that the average level density of Hamiltonians from the CUE, GOE and 
GSE ensembles take the form of a semi-circle [Meh04]. In other words, it says that 
on average the p.d.f. of the energy values of these systems is a semi-circle. 

As shown in the previous section Wigner’s Semi-circle law is in fact a statement 
about many-body hamiltonians, specihcally those where the order of the interaction 
k is taken to be equal to the number of particles in the system, m. Assuming that all 
particles interact simultaneously under the force of a potential was an assumption made 
implicitly in the canonical form of random matrix theory, because it led to important 
simplihcations in a challenging and very technical new held of mathematical endeavor. 
There was a second rationale as well. Measurements of particular nuclear energy levels 
in the 1950s and 1960s [RDRH60,DRHR60,DRHR64,HR51,MHR53] involved nuclear 
interactions with a high degree of random mixing. This was caused by high excitation 
energies due to the particles being accelerated within a synchrocyclotron. The 
interactions appeared to be so complex that a “black box” approach to studying the 
nuclei became a practical and reasonable way forward as it involved assuming the bare 
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(a) N=10 



(c) N=50 



(d) N=200 


Eigure 3.1; Plots of the eigenvalue distributions for N x N matrices sampled from 
the CUE ensemble with N increasing in each simulation. The entries of the matrix 
are taken to be complex gaussian random variables with variance 1. The symmetry 
condition placed on the matrices is H = W. 
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minimum of facts about the specific details of the potential H, beyond its symmetry. 
The most natural assumption was therefore to represent the physical hamiltonian as 
a hermitian matrix with some distribution (say a gaussian) determining the matrix 
elements; this was the black box. So apart from the mathematical simplihcation the 
initial approach to studying nuclear resonances, which resulted in the fonndation 
of canonical random matrix theory, was motivated by experiment. The fact that it 
implicitly implied k = m did not detract from the possibility that nature did indeed 
appear to behave this way, statistically speaking, in some complex systems. 
However, it soon became apparent that many-body potentials taking values of A; < m 
would need to be investigated as well. Indeed since the entire range 0 < k < m could 
plausibly model some experiment it became an imperative to hnd the statistics of 
the energy values for all of these possibilities, as this would be the only way to check 
whether or not the resnltant models predicted behavionr which is also measnred in 
actual experiments. This unihed approach to random matrix theory was proposed by 
K. K. Mon and J. B. French [MF75] and the nnihed ensembles determined by the 
modihed hamiltonian (2.1.86) were called the Embedded Ensembles. Hence, theoretical 
physicists and mathematicians hrst solved the case k = m and only later tnrned to 
the general case (still unsolved in its entirety) determined by allowing k to take any 
valne within the permitted range 0 < k < m. The next sections will explain how 
to calcnlate Wigners semi-circle law in two ways. One method will nse a complex 
technique called supersymmetry involving the liberal application of anti-commnting 
variables and their properties, which will be discussed beforehand. The second way 
will nse a simple method involving some basic diagrams and combinatorics. Both 
methods will be used to derive the semi-circle law of canonical random matrix theory, 
i.e. for the special case k = m. An attempt to calculate Wigner’s law for the general 
case k < m will be presented in later sections and will rely on new diagrammatic 
methods which are the primary topic of this thesis. There it will be seen that the 
analogue of Wigner’s Semi-Circle Law in the nnihed phase space of embedded RMT 
is, for many values of k, not even a semi-circle. 
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3.2 Diagrams 


In this section Wigner’s semi-circle law will be shown only for the case k = m with 


I —)■ oo and just one restriction on the symmetry of the potential, namely H = W. 
This is the canonical class of RMT known as the CUE ensemble (see section 2.1.2). 
For the purpose of the proof the quantum states containing m particles are assumed 
to consist entirely of fermions and these m-body states will be denoted by the greek 
symbols |/i), |z/), |p), |cr) et cetera. Each of the m particles in these states will be 
one of I single-particle occupation levels, with m I and the caveat that no two 
single-particle states in the same m-body state can be the same because the states 
are fermionic. Hence there is a total of A^ = orthogonal m-body states in the 
basis of the system. 


Notation 

The single-particle creation and annihilation operators are as before written as 


and Qj respectively with j = 1,..., /. With the intention of simplifying the notation 
dehne the shorthand expression 


3 {.jly ■ ■ ■ y jk) 


(3.2.1) 


and similarly for i. This furnishes a suitable abbreviation for the creation and 
annihilation operators which will now be written as 


~ • • • ®ji • 


(3.2.2) 


A useful corollary of this is the equality 



(3.2.3) 


Each of the ( M states in the basis can be written in the form 
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with |0) denoting the vacnnm state and the restriction 1 < ji < j 2 < • • • < jm < ^ 
since these are m-body states containing only fermions. From this it can be seen that 
for the special case k = m each m-body state can in fact be rewritten as 

at|0) (3.2.5) 

for some particular set j. The /c-body potential as derived in section 2.1.3 is written 
as 

Hk = ^ ■ ■ ■ o-h (3.2.6) 

with the stricter sum now running over fc-tuples with non-repeating states. With the 
notational abbreviations defined above this can now be rewritten in the following way 

Hk = Y^Vjia^jai. (3.2.7) 

ji 

Because present attention will be focused on the case where the potential is a hermitian 
matrix from the GUE ensemble, the only symmetry condition which Hk needs to 
satisfy is given by 


{fi\Hk\iy) = (z/|i7fe|/i)* (3.2.8) 

for all fi,u. This implies that 

Y,Vji{n\a\ai\u) = Y,v*j{u\a\ai\fi). (3.2.9) 

ji ji 

Hence for the specihc case k = m the expression (3.2.5) can be used to rewrite (3.2.9) 
as 

^nji(0|a^ataia|)|0) = '^v*j{0\a„a]aial\0). 
ji ji 


(3.2.10) 
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This grants a further simplihcation, since 

~ ^pq 

implying that (3.2.9) is equivalent to 


so that ultimately 




(3.2.11) 


(3.2.12) 


{n\H\u)* = = {i3\H\^). (3.2.13) 

It will become apparent shortly why this equality is going to be useful in subsequent 
calculations. 


3.2.1 Moments and Ensemble Averages 


The exact form of the average p.d.f. for the eigenvalues of the CUE ensemble will be 
calculated by finding the moments of the density function, the 2n-th moment being 
given by the expression 


/d2n 




(3.2.14) 


Fortunately (the moments would otherwise be inhnite), in the denominator there is 
always the normalisation expression 


1 

N 




1 _ 




2 


(3.2.15) 


with implicit summation over the repeated indices p, v and the overbar denoting 
the ensemble average (which is the value weighted by the probability that it occurs 
within the CUE ensemble). Yet to be determined is the probability distribution for 
those elements of the matrix which are not determined by symmetry. Following the 
norms of canonical RMT the gaussian distribution will be used, though practitioners 
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in the field do also study ensembles where this is not the case (e.g. [Bee94,IlW95]). 
The convergence of Random Matrix results for differing distributions of the matrix 
elements is referred to in the literature as universality [Meh04] . Assuming that all 
Vji are i.i.d. complex random variables sampled from a gaussian with mean zero and 
standard deviation 1 yields 

= N. (3.2.16) 

This permits another neat simplifcation; for the current case with k = m the moments 
are in fact given by 

To discover what the value of this trace is, re-write it explicitly in the form 


tr(hf^ ) — . . . Vj(2n)^{2n) (/i | U • (1) Cl j{l) Cl • (2) • • • Cl •( 2 n,) Clj{2n) |/i) (3.2.18) 


with summation over the repeated indices again implicit in the notation. Since these 
are complex gaussian random variables with mean zero the ensemble average denoted 
by the overbar will only be non-zero when it is made over even powers of the n’s. 
The ensemble average is the standard deviation which, to simplify notation, has 
been predetermined to be of unit value 


X;2 = 



1 . 


(3.2.19) 


The general case can be shown, given unit standard deviation, to give the ensemble 
average of every even power of v as 





= (2n - 


1 )!! 


(3.2.20) 


This, coincidentally, is the number of ways of uniquely pairing a set of 2n objects, 
which is also the sequence determining the moments of a gaussian distribution. It will 
be shown that the moments in this case are not in fact the moments of a gaussian 
distribution. 
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Products of Pairs of Random Variables 

Using just (3.2.20) equation (3.2.18) can be written as an expression determined by 
summing over only the product of the average of all possible pairs of the v. Specihcally, 
let {a} denote the set of (2n — 1)!! unique pairings of the sequence 1, 2,..., 2n with 
(t{x) taking the value of the integer paired with x for a given permutation a from 
the set {a}. Note that there is no x for which ct(x) = x. Then the average of the 
product of u’s from (3.2.18) can be rewritten in terms of averages of pairs of v 

2n 

'^j(l)i(l)‘^j(2)i(2) ■ ■ ■ Vj(2n)^(2n.) = ^ V j(x) l(x)Vjcy{x) la(x) . (3.2.21) 

CT X=l 

This result is used frequently in statistical mechanics and is sometimes referred to as 
Wick’s Theorem (even when it doesn’t involve creation and annihilation operators). 
Hence (3.2.18) becomes 




2n 

^ Vj{x)l(x)VjtT(x)j^a-{x) 

. cr x=l 


(/i|ci .(l)Cl^(l)Cl .(2)Clj(2) . . . O .(2„) UjCZn) l/i) . 


(3.2.22) 


Of note is that since k = m, operating on a state |/i) with annihilation operator Oi 
will give 

a*|/i) = (3.2.23) 

This, and the fact that the state |p) occurs on both the right and left hand side of the 
sequence of creation and annihilation operators of (3.2.22) means that this sequence 
must follow some path in state-space, hrst removing |/i) to give the vacuum state |0) 
then replacing it with some other state \t 2 n), removing this to give |0), replacing it 
with some other state, |t 2 n.-i) and so on, with the hnal step being a replacement of 
a vacuum state with the initial state |/r). This process of annihilation and creation 
occurs 4n times, with 2n annihilation operations and 2n creation operations. The 
result is the formation a path in state-space beginning and ending with the state |/r) 
which never removes a state which is not already there. Any terms in the sum of 
(3.2.22) which do not contain such a sequence will have a value of zero. Explicitly, 
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the rightmost component of the product given in (3.2.22) either takes the value zero 
or can be written in the form 

(hl«](i)ai(i)«](2)a*(2) • • • aj.(2„)aj(2n) |/i) = {^\alat^al^... = 1. (3.2.24) 

Hence 

(h'l®j(l)®i(l)®j{2)®i(2) • • • Q'j(2n)®i(2»^) IA'') = JT *^**^*’^4+1 (3.2.25) 

a b 

where H is identihed with /i and t 2 n+i is identihed with ti. The above allows (3.2.22) 
to attain the hnal form 

2n 

= E n (3.2.26) 

a a=l 

Each component of the product can be expressed as 

Dmtm+lD„t„+l (3.2.27) 

for some m, n. For uncorrelated and Vt^t„+^ the ensemble average will yield 

null, since the gaussian determining the p.d.f. of these elements is an even function. 
However, for correlated Vt^trr,+i and Vt„tn+i (3.2.27) will yield unity if and only if 

Vtmtm+l = (3.2.28) 

because the ensemble average will then yield the second moment of v which has been 
hxed to take the value 1. By (3.2.13) this is equivalent to saying 

+ (3.2.29) 


SO that (3.2.26) becomes 


2n 




cr a=l 


{a,) + l"*a + lto-(a) ' 


(3.2.30) 
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It will be shown in snbseqnent sections how this can be represented by diagrams. 
These will reveal that the product of delta functions giving (3.2.30) leads to a sequence 
of paths on the labels a = 1, 2,..., 2n and each path will be identihed with a partition 
of the set 1,2,..., 2n. Each partition will in turn be identihed with a diagram, and 
these diagrams will be counted to give (3.2.30) which can then be inserted back into 
(3.2.17) to give the moments of the level density for the CUE. 

3.2.2 Partitions 

Equation (3.2.30), and through it the moments of the level density (3.2.17), are 
related to certain combinatorial diagrams. This will be shown by observing that each 
instance of the sum of (3.2.30) can be mapped to a partition on the set of integers 
{1,2,..., 2n}. Using diagrammatic methods for the k = m case of canonical random 
matrix theory is not new [Kre72]. It will be seen that the moments of the level 
density are in fact given by the Catalan Numbers^, which occur in numerous other 
combinatorial problems outside of the held of random matrix theory. Eor these - 
some of which are much older than random matrix theory itself and have their own 
distinct literature - diagrams are frequently used. For 2?7, = 4 (3.2.30) gives 

_ 4 

'2^042^42*3^*3*4'^*40 n ^iaio-(a) + l^*ti+licr(a) ■ (3.2.31) 

(7 a=l 

the right hand side of which can be read as the “product of the average of all unique 
pairings of u”. Hence 

^tlt2^i3^4 ^^2^3^^40 
^tit2^i4ii ^i2i3^i3^4 

= Stlt3 + ^t2t4, + ^tlt2^t2t3^t3t4,^Uti (3.2.32) 


^the sequence defined by Cn = ;^(^") 
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Partitions as Cycles 

Identifying ti with 1, t 2 with 2 and so on, the components of (3.2.32) can be represented 
as cycles in the following way 

1.. .3 = (13)(2)(4) (3.2.33) 

3.3., = (24)(1)(3) (3.2.34) 

( 3 - 2 - 35 ) 


Note that where c* denotes the number of orbits in a given cycle with length i that 
for each cycle of (3.2.33 - 3.2.35) it holds that 


^^c^ = 4 = 2n. (3.2.36) 

i 

Each orbit represents a degree of freedom in the variables ti,t 2 , ■ ■ ■ (etc) of (3.2.31). 
The value of tr(if^) is given by implicitly summing over all values that each of these 
free variables can take, the result being that each orbit contributes a factor to 
the final expression. Hence the cycle (13) (2)(4) of (3.2.33) is to be identified with 

/ j \ 3 

ij and likewise for the cycle (24)(1)(3). The cycle (1234) however, consists only of 
a single orbit, so that it is identified with the value Hence 

tr(J?l) = 2(M +(') (3.2.37) 

\m J \m J 


so that the fourth moment for the GUE with k = m given by (3.2.17) is 


a _ to(i7^) 

jy2+l 



= 2 + 0 



(3.2.38) 


For 2n = 6 (3.2.30) reads 


^142^2*3 ^ 3 * 4^445 DstgDetl 
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+ Vtit2'Vt2t3 

T Clt2’^i2i3 

T Cit2^i3i4 ^t2i3^i4i5 ^*5^6^^6*1 
T '^tit2^i3i4 ^i2i3^i5i6 
T ^tlt2^i3i4 '^t2t3'Vt3t\ 

T 'i^tit2^i4i5 ^^2^3^*344 

+ Vtlt2'^t4t5 C2i3'^t5t6 
+ 'Vtlt2'^tit3 '^t2t3'^tetl Vt^tiVt^te 
+ Clt2'^i5i6 '^t2t3'^t3t4 Vt^ts'^tatl 

T "^42^3^i4i5 ^*3^4^46*1 

T 'Vtit2'^t3t(i '^t2t3'Vt3t\ '^1314^14,13 
T 1 ^tit2'^tQti '^t2t3'^t3t4 '^t4t3'Vt3tQ 

T ^tit2^i6ii ^^2^3^*445 

T Clt2'^i6il "^*2*3^*546 '^t3t4'^t4t3- 

Again identifying C with 1, t2 with 2 (etc), the components of (3.2.39) 
represented as the cycles 


^tit3^t3t3^t3ti 

113 A316 ^t415 A5 tg A4i 1 

Alta Asti A 4 i 6 

Ait 4 A2t3 Astg Asti Agti 
Aiti A2t3 A2t6 Agtg Agig 

Aiti A2t3 Ait2 Agtg Aitg 
Aitg A 2 t 4 

Aitg A2t4 A2tg A3t5 Aitg Aitg 
Alts A2t4Ait2 Agtg Alts 


(135)(2)(4)(6) 
(13654) (2) 

(13) (46)(2)(5) 
(14325) (6) 

(14) (2356) 
(14632) (5) 

(15) (24)(3)(6) 
(153)(246) 
(1542)(36) 


(3.2.39) 
can be 
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= (16524) (3) 

^tits^t2ti^t3t4^tlt3^t4,te = (1643) (25) 

<^4245^41*2^*3*6^^3*5 ~ (16352) (4) 

^*2*6^*2*4^*4*6 ~ (246)(1)(3)(5) 

^*2*6 ^*2*5 <^*3*4 <^*3*6 <^*4*5 ~ (26345) (1) 

^ (26)(35)(1)(4). (3.2.40) 


This time the cycles obey 


y^ici = 6 = 2n 


(3.2.41) 


and since the valne of each term is the nnmber of orbits in the cycle times 
gains 



one 



(3.2.42) 


so 



(3.2.43) 


For 2n = 8 (3.2.30) reads 


tr (FT^ ) '***1*2^2*3^344'***4*5'**46*7^7*3'44*8*1 • 


(3.2.44) 


Now there are (2n — 1)!! = 105 ways to pair the n’s, each giving rise to a unique 
partition on the integers 1, 2,..., 8. This time there are just 14 cycles determining 
the leading behaviour of the 8-th moment a,s N ^ 00 , given by 


(2468)(1)(3)(5)(7) 

(1367)(2)(4)(6)(8) 


(3.2.45) 


(15)(24)(68)(3)(7) 

(17)(26)(35)(4)(8) 

(28)(37)(46)(1)(5) 
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(13)(48)(57)(2)(6) (3.2.46) 

(135)(68)(2)(4)(7) 

(248)(57)(1)(3)(6) 

(137)(46)(2)(5)(8) 

(246)(17)(3)(5)(8) 

(357)(28)(1)(4)(6) 

(268)(35)(1)(4)(7) 

(468)(13)(2)(5)(7) 

(157)(24)(3)(6)(8). (3.2.47) 

Note that the cycles (1357)(2)(4)(6)(8) and (2468)(1)(3)(5)(7) of (3.2.45) are equiv¬ 
alent up to a constant mod 9. In other words the orbit (1357) with each element 
incremented by 1 mod 9 is identical to the orbit (2468) and similarly the orbits 
(2) (4) (6) (8) incremented by 1 mod 9 are equal to (3)(5)(7)(1). The same applies to 
the four cycles of (3.2.46), each being equivalent to eachother up to an increment of 1 
mod 9, 2 mod 9 or 3 mod 9. Finally, the eight cycles given by (3.2.47) are equivalent 
up to an increment of i mod 9 where 1 < i < 7. Assigning each orbit of each cycle 
the value as before and dividing by the normalisation constant A^”+^ gives the 
eighth moment for the CUE as 

A = ^ = 14 + O (1) . (3.2.48) 

The process for calculating the moments could continue indefinitely in this way. But 
what is the underlying pattern? And in what form can this pattern be codihed in 
order to hnd a general expression for all moments of the k = m scenario? In the 
next section the cycles introduced above will be represented as diagrams. This will 
illustrate the possibility of transferring the problem (of calculating the moments) 
into a “diagram space” and using the new space to make calculations. The method 
will act as a precursor to a much more complex diagrammatic technique which will 
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be used in subsequent chapters to calculate moments for random matrix ensembles 
where k <m. 

3.2.3 Polygons 

The first method for representing (3.2.30) will be to simply draw the cycles from the 
previous section as segments of 2n-sided polygons. It will be seen that the leading 
order behaviour of the moments, being the cycles with n + 1 orbits as shown in 
section 3.2.2, are given by those partitioned polygons where (i) no partitions intersect, 
(ii) no two isolated vertices neighbour eachother and (iii) any two vertices which 
can be connected by a non-intersecting edge are connected as such. For the 4-th 
moment there are just two cycles with three orbits, namely = (13)(2)(4) and 
5t2u = (24) (1) (3) which are illustrated in Fig 3.2. 



Figure 3.2: Plot of polygons for fourth moments with the left hand side giving an 
illustration of the cycle (13) (2) (4) on a 2n = 4 sided polygon and the right hand side 
illustrating the only unique rotation of this cycle, which is (24) (1) (3). 




Figure 3.3; Illustration of the cycle (135)(2)(4)(6) and it’s unique rotation 
(246) (1) (3) (5) given by incrementing each element by 1 mod 7. 
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Figure 3.4: Illustration of the cycles (13)(46)(2)(5) (left), 
(26)(35)(1)(4) (right) of (3.2.40). 




Figure 3.5: The cycles (2468) (1) (3) (5) (7) (left) and (1357) (2) (4) (5) (8). 
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Figure 3.6; The four equivalent cycles (up to addition modulo 9) given by 
(24)(15)(68)(3)(7) (top left), (17)(26)(35)(4)(8) (top right), (28)(37)(46)(1)(5) (bottom 
left) and (13)(48)(57)(2)(6) (bottom right). 
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6 7 



3 2 

(a) (135)(68)(2)(4)(7) 



(d) (246)(17)(3)(5)(8) 
6 7 



3 2 

(g) (468)(13)(2)(5)(7) 


6 7 6 7 






6 7 



Eigure 3.7: Illustration of the eight cycles of equation (3.2.47) each with the same 
shape, but rigidly rotated about the origin. 
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For the sixth moment the contributing cycles from (3.2.40) are illustrated in Figs 
3.3 and 3.4. Finally the 14 contributing cycles for the 8-th moment are illustrated 
in Figs 3.5 - 3.7. These illustrations prompt the following conclusion; the 2n-th 
moment of the level density of the GUE for the canonical case k = m is equal to the 
number of was of correctly partitioning a 2?7,-sided poygon, where the conditions on 
the partitions are (i) no partitions intersect, (ii) no two isolated vertices neighbour 
eachother and (iii) any two vertices which can be connected by a non-intersecting 
edge are connected as such. It turns out this is an old problem and has been studied 
before by the mathematician Germain Kreweras in the early 1970’s [Kre72]. In the 
next section the problem of counting the number of correct partitians of a 2?7,-sided 
polygon will be translated into another old problem; Dyck Words. In this way it will 
be shown that the 2n-th moment of the GUE with k = mis given by the sequence of 
Gatalan Numbers 1, 2, 5,14,42,132,429,1430,4862,16796,... defined by 


/^2n 


1 

n -I- 1 



(3.2.49) 


3.2.4 Dyck Words 

Each instance of the sum (3.2.26) gives a cycle and each cycle can be associated with 
the partitioning of a 2n-sided polygon. These in turn can be translated into Dyck 
Words. Equations (3.2.33, 3.2.34) give the two cycles of the fourth moment as 

(13)(2)(4) 

(24)(1)(3) (3.2.50) 


These are translated as follows 

(i) Write the sequence 123..2n in order without brackets 

(ii) A number with a left bracket in the partition notation is also given a left bracket 

(iii) A number with a right bracket in the partition notation is also given a right 
bracket 
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(iv) Remove all numbers, leaving only sequences of opening and closing brackets 

(v) Replace each opening bracket with a X and each closing bracket with a Y. 

The result is a Dyck Word. As an example take the cycle (13) (2) (4). For this 2n = 4 
cycle the ordered sequence 1234 hrstly becomes (1(2)3) (4) because the partition 
(13)(2)(4) has a left bracket for the integers 1,2 and 4 and a right bracket for the 
integers 2,3 and 4, as does the ordered sequence of numbers and brackets (1(2)3)(4). 
This in turn becomes (())() by removing all numbers, and hnally by (iv) the result is 
a Dyck word of length 2?7, -|- 2 = 6 

{1?,){2){A)=XXYYXY (3.2.51) 

by replacing each ( with X and each ) with Y. Similarly the cycle (24) (1) (3) becomes 
(1)(2(3)4) becomes ()(()) giving the 6 letter Dyck word 

{2A){l){?,)=XYXXYY. (3.2.52) 

It should be clear that this translation can be performed for every cycle (3.2.26) 
determining the moments of the level density. However the converse is not true. Not 
every Dyck word of length 2n refers to a permitted cycle of length 2n -|- 2. The Dyck 
words of length 6 for example are given by 

XXXYYY, XXYXYY, XYXYXY, XXYYXY, XYXXYY (3.2.53) 

whereas the fourth moment of the level density results in the cycles (13) (2) (4) and 
(24)(1)(3) (see 3.2.33, 3.2.34) corresponding only to the last two Dyck words of (3.2.53). 
Hence the amount of information contained in a permitted partition of 2n elements 
is not equal to the information represented by a 2n -|- 2 Dyck word. This can be 
appreciated by observing that the cycles (3.2.33- 3.2.34), (3.2.40) and (3.2.45-3.2.47) 
with polygon representations given in their respective diagrams Figs 3.2 - 3.7 are all 
completely determined by the cycles containing all the odd (or equivalently even) 
integers of the sequence. In other words, given the restrictions of section (3.2.3) 
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each diagram can be reproduced knowing only the cycles involving odd integers, or 
conversely knowing only the cycles with even integers. For the fourth moment the 
relevant cycles (3.2.33- 3.2.34) become 

(13)(2)(4) = (13) 

(24)(1)(3) = (1)(3). (3.2.54) 

In this way each polygon for the 2n-th moment is associated to a partition of n odd 
numbers as illustrated in (3.2.54). Because these sequences consist of odd numbers it 
helps to make the identification 1—)-l,3—?-2,5—>-3, ...,2n — 1—^-nso that (3.2.54) 
ultimately is mapped to 


(13)(2)(4) = (13) = (12) 

(24)(1)(3) = (1)(3) = (1)(2). (3.2.55) 

where the information given by the partitions on the l.h.s of the relation is equivalent 
to the information contained on the r.h.s and both can be used to reproduce a 
contributing partition for a 2n-sided polygon. The cycles for the sixth moment can 
similarly be rewritten as 


(135)(2)(4)(6) = (135) = (123) 
(13)(46)(2)(5) = (13)(5) = (12)(3) 
(15)(24)(3)(6) = (15)(3) = (13)(2) 
(246)(1)(3)(5) = (1)(3)(5) = (1)(2)(3) 
(26)(35)(1)(4) = (35)(1) = (23)(1). 


(3.2.56) 


It can be seen that each of the permitted cycles contributing to the leading order 
term of the 2?7,-th moment is equivalently represented as a non-crossing partition of n 
integers, so that the number of cycles which give the value of the 2n-th moment is 
equal to the number of non-crossing partitions of the set {1,2,..., n\. This implies 
that summation is identical to the summation of Dyck words of length 2n, not 2n + 2. 
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This is a very famous, beautiful and old result which has long been solved. As 
illustrated above, it has numerous forms including the ennmeration of non-crossing 
partitions of {1, 2,..., n}, as well as the partitioning of polygons in varions ways. 


Catalan Numbers 


Every Dyck word of length 2n can be represented as monotonic paths drawn in a 
n X n grid. The number of Dyck words of length 2n is given by the n’th Catalan 
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Fignre 3.8: Plot of the Dyck Paths. 


nnmber 


C'„,= 


n -|- 1 


(3.2.57) 


which is equivalent to saying that the level density of (4.2.5) is a semi-circle, since 
these are the moments of a semi-circular probability density function. In other words 


2 

vrr^ 



dx = Cn- 


(3.2.58) 


We will briefly look at one way of showing that the nnmber of Dyck words of length 
2n is indeed given by 3.2.57. There are many cute proofs to this and the one presented 






























3.2. DIAGRAMS 


51 


next is to illustrate the fact that a diagrammatic method can swiftly yield the semi¬ 
circle law. 


Proof 

The proof begins with some monotonic sequence involving n vertical steps and n 
horizontal steps which do indeed cross the horizontal line, and therefore do not 
constitute an acceptable Dyck word. Identifying the hrst node lying above the 
diagonal, and drawing a line through this point which runs parrallel to the diagonal 
gives the picture illustrated in hg 3.9a where all steps taken after this point are 
coloured in red. Dehning the number of horizontal steps taken up to this point as fc, 
it follows that the number of vertical steps taken is fc -|- 1, since there is one arrow 
peaking up above the diagonal. Now reflect the remaining n — k horizontal steps. 
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(a) A non-Dyck word on a n X n grid. 
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(b) Monotonic path on {n — l)x{n + 1) grid. 


Figure 3.9: Illustration of the transformation showing that every non-Dyck word on 
a n X n grid can be mapped to a monotonic path on a (n — 1) x (n -|- 1) grid. 


and the remaining n — k — 1 vertical steps in the line drawn parrallel to the diagonal 
passing through the hrst node above it. This yields the diagram shown in hg 3.9b. 
The hipping step has turned n — k horizontal steps into vertical steps and n — k — 1 
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vertical steps into horizontal steps. Hence the resulting picture is a monotonic path on 
a {n — k — l + k) x {n — k + k + 1) grid. In other words, any monotonic path from the 
initial n x n grid which is not a Dyck word is a monotonic path on a (n — 1) x (n + 1) 
grid. The number of Dyck words of length 2n is then given by 


'2n\ / 2u \ _ 1 /2n\ 
n j \n + ij n + l\n J 


(3.2.59) 


since we hrst hnd all ways of making a monotonic path from n vertical steps and 
n horizontal steps, and then remove all monotonic paths making (n — 1) horizontal 
steps and (u -|- 1) vertical steps. □ 

With relative simplicity it has been shown how to calculate a random matrix 
theory result for the k = m case. But what about when k < m? In later chapters it 
will be shown how diagrams can once again be utilised in calculations of the moments. 
Whereas for k < m these diagrams will appear and function differently to what has 
been shown above, for A; = m they become equivalent to the representation of cycles 
on polygons as shown in section 3.2.3. 


3.3 Supersymmetry 

In this section it will be seen how an existing technique involving grassmannian 
variables and an innovation by [BRWOla] can be used to show that the level density 
of the embedded (a.k.a unified) CUE ensemble is a semi-circle, as may be expected. 
However, this result will only cover the domain m <2k with I —)■ cx). In the context 
of unified RMT the result covers the volume of the phase space {k > y,m,cx)}. In 
the next chapter a new method will be introduced to calculate statistics of the unihed 
RMT phase space for values of 2A: < m as well. It will be seen that in this portion 
of the volume of the phase space the analogue of Wigner’s semi-circle law is not a 
semi-circle! 
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3.3.1 Anti-commutative Grassmann variables 

The proof of Wigner’s semi-circle law for fc-body potentials begins nsing the same 
syperanalysis as nsed for the m-body case. A good acconnt of this can be fonnd 
in [HaalO, AMPZJ94]. One begins with dehnitions which describe the form and 
behaviour of anti-commutative variables, which are non-numerical “things^’ with the 
following algebraic rules attached 

P: Q = ^ ■ ■ ■ d (3.3.1) 

which implies, among other things, that = 0. Dehning := to be the 
complex conjugate of ^p, where is another grassmannian variable which therefore 
must also obey (3.3.1), and additionally enforcing the rules 


(&«.)• = C;C,C (3.3.2) 

and 

Cp* = -^p (3.3.3) 

it follows that {^p^p)* = ^p*^p = ^p^p- In other words, what has been called “complex 
conjugation” leaves the grassmannian ^p^p unchanged. In this sense it can be thought 
of similarly to the length of a standard complex variable. Finding a complex analogy 
for every grassmanian property, however, is not the goal. They are ultimately very 
different, distinct objects. Another property of grassmannians which illustrates this 
well is that they do not have an inverse. Namely, assuming that exists gives 
^P^^P^P^ = -^p\^p^Q = -^p^ and likewise = (^p ^^p)^p ^ = ^p ^ so that 

^p^ = —^p^- Namely, it is not possible to abide by the original assumptions, (3.3.1), 
(3.3.2) and (3.3.3) without also contradicting the existence of an inverse. Analogously 
to the way in which a real or complex C“ function can be expressed in the form of a 
Taylor series f{x) = a+ bx + cx^ + ... any function of purely grassmannian variables 
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can be expressed in the form 


E (3.3.4) 

since = 0 for any n >2. Having defined the form of grassmanian functions it is 
now possible to describe a formal meaning for differentiation and integration. That 
is, since grassmannian objects are different to numerical objects it is necessary to 
redehne what it means to differentiate and integrate. Moreover, there is no constraint 
to make these dehnitions in any way analogous to the numerical case, but rather the 
dehnitions are only required to be commensurate with (3.3.1), (3.3.2) and (3.3.3). 
The dehnitions are 


/ = 1 

(3.3.5) 

1 = 0 

(3.3.6) 


where the differential d^i obeys the same commutation relations as a normal grass¬ 
mannian variable. Grassmannian differentiation is dehned analogously to real dif¬ 
ferentiation, namely d/d^i{a^i) = a for any commuting variable a. Moreover, the 
operator d/d^i obeys the same commutation relations as grassmannian variables so 
that 

d d d d 

dii dij d^j dii 

which yields jdQ = 0 and dld^i{^j) = —{ij)d/d^i for i ^ j. It follows immediately 
that 

= (3.3.8) 

and likewise 

J dU...d^,F = f{l,l,...,l). (3.3.9) 

That is, added to the already remarkable properties of anti-commuting variables is 
the fact that the operation of differentiation is identical to that of integration. 


(3.3.7) 
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3.3.2 Determinants 

From this comes 


'1...7V 


1...N 


D = I [E exp{- 

j / ik 

I...N \ 

n dCjdQ exp{-CAC) 


3 


(- 1 ) 


N 


n n’Ai I 


3 


^3 ‘’J/ JY! 
1 


n HA] I 


(3.3.10) 


since by (3.3.6) the Y’th term is the only one contribnting a non-zero term to the 
integral. There are (Y!)^ terms in the prodnct of the snm. N\ for the nnmber of 
ways a non-zero term can be chosen, times N\ for the number of different orderings 
for such a term. Hence 

‘^ = 1 ('if4 (3.3.11) 

where i and j are permutations on the indices. Since each permutation on i gives 
the same contribution it follows that 


T...Af 


^=y n dipc^ E a,p.a,p, ... a 

= AiPiAi^p^ ... Ai^p^ 

p 

= det A. 


ijvPjv 


(3.3.12) 


So using grassmannian variables the determinant of a matrix can be expressed as 
the integral of an exponential function of anti-commuting variables. This will prove 
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important later on. Another important is the observation that 



(3.3.13) 


where d'^z = dRe{z}dIm{z} and it has been assumed that A is hermitian, and 
making the variable transformation Qk = {Uz)k for A = t/lAC. The determinant of 
the Jacobian relating the differentials d‘^[Z] and d?‘[Q] is subsequently unity. Hence 

/ ^ (n -p E = nW “ so- 


3.3.3 The Generating Function 

It will now be shown how the Green’s function can be expressed in terms of a 
generating function. The grassmannian variables and the resulting expressions for 
determinants found in the preceding section will then be used to express the generating 
functions in terms of gaussian integrals. This process of changing the form of the 
determinant into a gaussian integral using functions of grassmannian and complex 
variables is called the method of supersymmetry. The level density of the system is 
p{E) = ^Im G{E~) where the Green’s function G{E) := j^Tr Introducing 

the generating function 


Z{E,j) 


det {E-H) _ det i{E - H) 

det {E-H -j) ~ det i{E - H - j) 


(3.3.15) 


and observing that 


dj 


det {E-H- j) 


Ji=o 


= ^exp [-Tr{ln{E -H - i))]j=o 


d 

= —exp 
dj 


-Tr{ln{E - H)) + ]Tr{E - H)-^ + 0{]^) 


3=0 
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1 

det (E - H) 


Tr 


- 1 ■ 
.E-H. 


gives the Green’s function as 


(3.3.16) 


G{E) 




Z{E,j) 



(3.3.17) 


Given that p{Ei, E 2 ) = p{Ei)p{E 2 ) it is now possible to write the two point generating 
function as 


Z{E^,E+,E^,Et) 


detiE^ - H)det{Et - H) 
det(Ef - H)dei{Et - H) 

^ dei[i{E^ - H)]dei[i{Et - H)] 

^ ’ det[i(Ef - H)]det[-i{Et - H)] 


(3.3.18) 


where using the notation of (3.3.15) gives E^ = (Es — j) — iei, E 2 = (E 4 — k) + ie 2 , 
E^ = E 3 — iei and E^ = 7^4 + ie 2 , implicitly with 64,2 —t 0 so that 


G{E^)G{Et) 


1 d^Z{E^,Et,E^,Et) 


Y2 dEidE2 


E\=E^,E2=E4, 


(3.3.19) 


Plugging (3.3.12) and (3.3.14) into the above immediately yields 

X exp[f(4(77 - E 7 )zi + 4(77 - E^)pi)] 

X exp[i(-4(77 - E:^)z 2 + 4(77 - E^)p2)] 


(3.3.20) 
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where Z\^Z 2 are A^-component complex variables and rii,ri 2 are Wcomponent grass- 
mannians. Now for brevity define 





fl 
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[V2 J 


lo 

0 

0 



(3.3.21) 


so that d^*d^ = d[zi,ril,rii]d[z 2 ,ri 2 Tri 2 ], the matrix L is 4A^ x 4A^ with La '■ = 

Laa and the generating fnnction becomes 


Z = (-1)^ j d^*d^ exp ^ ^lLa[H - Ea]^c 
Hence, for the ensemble average the aim is to caicniate 


(3.3.22) 


IV := exp(zJ2^aLaff^c 


= / d[v] exp i ^ ^aLaH^a - 


2n2 

Jl O 


= J d[v] exp [I'^aT . (3.3.23) 

V ° ji 

Making the change of variables Vji —)■ vji — ivlJ2^a^La{p\a\aj\h')^au (notice the 
change in order of the creation and annihilation operators in the second term) and 
integrating over v gives W as 




(3.3.24) 


Using the method of expansion in the eigenvalues detailed in [BRWOla] one has 
A^upa ■= {l\VW){p\VW) = ^ Substituting this into (3.3.24) and 

taking caution with the anti-commutative properties of the grassmannian components 
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of $ grants the following equality 


fU = exp I Y. LaL'^Qa/^Qpc 


OLf3sa 


where 


L' = 


/ 1 0 0 0 \ 

0-10 0 
0 0-10 
0 0 0 -1 


and with A® := a/A*. 


(3.3.25) 


(3.3.26) 


3.3.4 Supertrace and Superdeterminant 

By writing out the 16 components of (3.3.25) explicitly it can be rewritten as 


W = exp 



with Q 


LaQajS and the supertrace Str defined as 


StrF = TrA — TrU 


(3.3.27) 


(3.3.28) 


where F 


A 

C 


B 

D 


is a supermatrix, the matrices A and B consisting only of 


commuting variables (such as complex numbers), and the matrices B and C consisting 
only of grassmanian variables. It is now possible to introduce an additional dummy 
variable by noting that up to a normalization factor 
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Then writing out the expression component-wise and rearranging it can be seen that 


Str aQ = ^^L{saf{vlX^C^'^h)^ 


(3.3.30) 


where s = diag{l 2 , —I 2 ) = so that up to a normalisation factor the ensemble 
average of the generating function Z is given by 


exp ( -yStr + ^^[iL{s(Tf {vlyC^^^h) - iLE]^ 


(3.3.31) 


to which the identity / exp(—'h^’Fd)) = (SdetF) ^ can now be applied, with 

the superdeterminant Sdet being dehned by 


Sdet = 


det(Zl - CA-^B) 
detA 


(3.3.32) 


Noting that this will be a AN x AN superdeterminant (denoted in bold to distinguish 
it from the 4x4 superdeterminant) and again ignoring normalisation constants 

(3.3.31) can now be written as 

J (Sdet[L(s(T)^(noA*C®“) — L£'])“^exp • (3.3.33) 

Using InSdet F = Strln F and factoring out L in addition to noticing that 


exp(StrlnL) = exp Str 


f 0 0 0 0 \ 

0 ZTT 0 0 

0 0 0 0 

0 0 0 0 


= exp(fA^7r) 


(3.3.34) 


and ignoring the ±1 implied by (3.3.34) it can be concluded that Z without normal¬ 
ization can be written as 


exp ( -yStr + Strln (e - 


(3.3.35) 
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Finally, dropping the transpose due to the diagonal nature of the term in the 

4x4 sense and reconciling the dimensions of the two supertraces via Strln = trStrln 
yields 

Z oc y exp Y[-Str + —trStr ln(crsUoA®C*“ — E)]^ . (3.3.36) 


3.3.5 Saddle-point Approximation 

A saddle-point approximation will now be applied over the a variables. Minimising 
the argument of the exponential yields 


a™ 



V X-E ) 


(3.3.37) 


where y = Multiplying across by and summing over sa as 

well as expressing the trace in component form gives 

t E (-E - A*C“C“ = vl E(S - xUWWjWi?)- (3.3.38) 

^ sapa per 


Assuming that y is of the form rl for some constant r and given that all the 
eigenvectors C^°' are traceless save for = 5^,^ gives, self-consistent with the prior 
assumptions on y, that 


X = 


u^A° 

E-x 


(3.3.39) 


where it is shown in [BRWOla] that 


A°(fc) 


mA /l — m + k\ 
k)[ k )■ 


Finally, solving this for y yields 


(3.3.40) 
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with t±{E) = i Y 1 — (^ 2v^xo '} ■ Pluggiiig this into (3.3.37) implies = t and 

cr®“ = 0 for all s > 1. Construct the supermatrix a as 

cr° = diag(r+(Ei),r_(E2),r+(E3),r_(E4)) (3.3.42) 


so that the supertrace for cr^ doesn’t yield zero and the advanced and retarded 
components of the generating function are matched. Then by setting E 2 = E^ 
in (3.3.36) the two-point generating function reduces to a level density generating 
function Z{Ei^E^) proportional to 


/ 


exp 






exp 


{-N [i,i(A",.y+(s.) 


£i)-lii(AVr+(B3)-£3)]). 

(3.3.43) 


Taking the derivative with respect to Ei as per (3.3.19) and selecting the imaginary 
component of the resulting Green’s function gives the embedded Wigner’s semi-circle 
law for the fc-body EGUE as 


lim p{E) oc 

N^OO 



(3.3.44) 


2 

Due to a technicality involving the identity ^ J2sa fh® saddle-point 

approximation converges to this result only for 2A; > m. In the next chapter a simpler 
method will be introduced, which avoids the grassmannian gymnastics illustrated 
above and can be used to find the moments of embedded RMT systems for all k. 
















Chapter 4 


Many-Body RMT 


4.1 RMT vs. Embedded RMT 

Many-body random matrix theory is the application of random matrix concepts to 
the stndy of the fc-body potentials (see section 2.1.3). This is a superset of random 
matrix theory with hamiltonian 

Hk = ^ (4-1.1) 

l<zi < ... 

and coefficients taking the value of a random variable. By determining the 

coefficients from a probability density function the static potential (4.1.1) generates 
an ensemble of matrices which can be studied statistically, just as in the case of 
canonical RMT. As shown in 2.1.4 the canonical form of RMT coincides with just 
a single point k = m in the phase space of the unihed random matrix theory, also 
known as embedded RMT. One of the purposes of unffication is to study all of these 
random matrix theories as a group, rather than creating new methods and theorems 
for each particular case. In this sense embedded RMT treats all randomised fc-body 
systems equally, allowing k to take any value in the range 1 < k < m. While there is 
a vast literature for the statistics and theorems of canonical random matrix ensembles 
{k = m) including Wigner’s semi-circle law described in chapter 3, little was known 
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about the properties of unified RMT ensembles [1 < k < m) including the analogue 
of Wigner’s Law for these embedded systems until relatively recently. In fact it turns 
out that for some values oi k < m the analogue of Wigner’s semi-circle law is not 
even semi-circular. 

While some progress has previously been made towards calculating various statis¬ 
tics for embedded ensembles [MF75, BRWOla], the moments for the level density 
for these systems when k < 'j were unknown beyond the fonrth moment nntil the 
method invented by the anthor was introdnced. The cnrrent chapter will introdnce 
this new method, called the method of particle diagrams, and illustrate how it can be 
nsed to calcnlate the fonrth, sixth and eighth moments of the level density for the 
eGUE ensemble. Recall from see section 2.1.3 that the eGUE ensemble represents the 
hamiltonian of non time-reversal invariant quantum systems of m particles interacting 
nnder the force of a fc-body potential. Interestingly, the moments calculated in 
this way point to a convergence in the statistical behaviours across a wide range of 
many-body hamiltonians of a similar form, albeit nsed independently to stndy the 
statistics of quantum spin chains, graphs and hypergraphs [HMH05, ES14, KLW14]. 

In the context of embedded RMT the three symmetry groups introduced by Dyson 
for canonical RMT become the embedded GUE (eGUE), embedded GOE (eGOE) 
and embedded GSE (eGSE). For each of these classes the hamiltonian (4.1.1) obeys 
the same symmetry rules as in the canonical case even though the resulting statistics 
can be completely distinct dne to the changing strnctnre of H itself. The only 
non-deterministic components of the potential are the random variables 
With the foresight that their symmetry properties will be reqnired in calculations 
later, the next sections will investigate the second moments of Vji for the eGUE, 
eGOE and eGSE. This will also illnstrate how the restrictions of symmetry from 
each of Dyson’s three groups affect the symmetry properties of an embedded fc-body 
hamiltonian. 
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4.1.1 Second Moments of the Random Variables 


R was shown in section 2.1.3 that the hermitian potential describing the energy 
of a system of m fermions interacting nnder the inflnence of a fc-body potential is 
Hk = j: <...in■ •* l^he SGCOTid TnoTn/STits of the random 

l<i^ <. ..<i^ </ 

variables = Vji are dehned as For the method proposed in this 

thesis for calculating the moments of the level density, the socalled second moments 
form the first of several essential ingredients. The normalised 2n-th moments of the 
level density are given by 




(4.1.2) 


for which the numerator can be re-written as 


) — Vj(l)j^(l)Vj(2)l{2) . . . 11^(271) j(2n) (/i I Cl .(2) flj(2) . . . Cl .(2„) Clj(2n) | /i) . (4.1.3) 


Even before seeing the forthcoming calculations for the moments of the eGUE it can 
be appreciated that in order to calculate (4.1.3) it is necessary to also calculate the 
ensemble average 

Vj(l)i(l)Vj{2)i(2) . . . Vj(2n)i(2n). (4.1.4) 

It has been noted before in section 3.2.1 that the second moment thpvF = 0 for 

J j 1 

uncorrelated Vji,Vj'^' whereas for correlated with Vji = nV/ the second 

moment is unity (by virtue of normalisation). This and the additional knowledge 
that in general for a gaussian random variable 

J = (2n — 1)!! (4.1.5) 

gives the ensemble average of the product of random variables (4.1.4) in terms of the 
product of the average of all possible pairings of the random variables 



Vj(l)i(l)Vj{2)i(2) . . . Vj(2n)i(2n) 


2n 

'y ) Vj(x)j^(x)Vjcr(x)j^cr(x). 

cr x=l 


(4.1.6) 
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This result is used frequently in statistical mechanics and is sometimes referred to as 
Wick’s Theorem (even when it doesn’t involve creation and annihilation operators). 
Since the hrst step for calculating the higher moments is to calculate these second 
moments the next sections will explain how to calculate the quantity VjWjW 
eCUE, eGOE and eGSE. 

4.1.2 Unitary Symmetry {(5 = 2) 

It was shown in section 2.1.2 that the CUE ensemble refers to the set of random 
matrices obeying hermitian symmetry 

Hk = Hl (4.1.7) 

which are the class of hamiltonians referring to time-reversal invariant fermionic 
qnantnm systems. Since by the Panli exclnsion principle many-body fermionic 
quantum states cannot contain repeated single-particle states (4.1.1) can be rewritten 
as the restricted sum 

Hk 'y ) (4.1.8) 

As seen before in section 3.2 there is a useful vectorised abbreviation of this given by 

Hk = Y^Vjia]ai (4.1.9) 

ji 

where j = {ji, ■ ■ ■ ,jk) (similarly for i) and the corrolary Oj = ... aj-^. In a fermionic 

state space components of the sum where the vectors j or i contain repeated indices 
will anyway be zero, so restricting the snmmation is not actnally necessary. To 
generate an ensemble of matrices from this potential each of the variables Vji in the 
sum is taken to be a complex valued random variable with a gaussian p.d.f. of mean 
zero and unit variance. The system described by the hamiltonian is assumed to have 
some number I > m non-degenerate single particle levels so that the dimension of the 
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space is N = A matrix symmetrised with the condition (4.1.7) obeys 

(/i|i7fc|z/) = (z/|i7fc|/i)* (4.1.10) 

for all qnantnm states |/i) and \i>). This implies that 

(4.1.11) 

ji ji 

and matching coefficients for each term in the sum gives 

v'j, = Vij. ( 4 . 1 . 12 ) 

The coefficients Vji not related by (4.1.12) are then uncorrelated i.i.d. gaussian random 
variables with mean zero and variance 1. This means that for uncorrelated vji and 
^j'i' 

VjiVj'i' = ^ / exp dvjidvj'^' = 0 (4.1.13) 

since the gaussians are even functions, making the total argument of the integral odd 
with respect to the integration variables. Likewise if Vji = v^i^' then Vji and are 
complex and their product yields two even functions of opposite sign, again giving a 
zero average. However in the hnal scenario with 

Vji = V*.,., (4.1.14) 


the second moment becomes 



nl^exp 



dvji = 


(4.1.15) 


from which it can be directly concluded that 






68 


CHAPTER 4. MANY-BODYRMT 


This is the same identity calculated in section 3.2 for the case k = m. It will be an 
important identity in subsequent sections for calculations made in the space of unified 
random matrix theory. The reason it will play such a fundamental role in calculations 
is based on the idea behind Wick’s Theorem; each term in certain summations of 
products will be expressed in terms of summations of products of pairs, for which 
(4.1.16) will become an important aid. Analogous identities for the eGOE and eGSE 
will be derived next. 


4.1.3 Orthogonal Symmetry (/3 = 1) 

The embedded gaussian orthogonal ensemble (eGOE) refers to the set of potentials of 
bosonic time-reversal invariant quantum systems, as discussed in section 2.1.2. The 
hamiltonian 

Hk = ^ (4.1.17) 

l<il<...<ik<l 

sums over m-body bosonic states which are permitted to have repeated single-particle 
indices. The number of bosonic m-particle states containing repeats of z < m unique 
single-particle states is dimension of the state space is then just the 

sum of all these 




(4.1.18) 


As explained in section 2.1.2 the symmetry properties imposed on Hk imply that it 
consists only of real values, while also obeying the hermitian symmetry satisfied by 
matrices in the eGUE (see previous section). Hence for {3 = 1 




(4.1.19) 


Vji = Vji = Vij 


Matching coefficients gives 


(4.1.20) 
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so that with Vji = or Vji = v*,., the value of vjjvy^ will be 1 while for uncorrelated 
Vji, Vji-' the value of will be zero. There will be no cancellations in the integral 

since these are real random variables, so the product of Vji with itself yields only a 
single term. There are no cross terms as there would be if Vji were complex valued as 
in the case of the eGUE. Hence for the eGOE 


(4.1.21) 



Analogously to the eGUE, this identity is important for calculations of the moments 
of the level density for taken from the eGOE because the resulting expressions 
involve summations of products of Vji which can be expressed as summations of 
products of pairs of Vji, for which (4.1.21) then plays a central role. 

4.1.4 Symplectic Symmetry (/3 = 4 ) 

For /3 = 4 the hamiltonian must satisfy both the hermitian symmetry of (4.1.7) and 
the symmetry implied by the fact that the potential takes the form of (2.1.33). To 
begin, dehne the pairs of entries in the upper(lower) diagonal plane of (2.1.35) related 
by complex conjugation as blue pairs, and those related by complex conjugation times 
— 1 as red pairs. The additional symmetry enforced by (2.1.33) implies that for every 
blue matrix element there is a map to the element with row number and column 
number a distance +1 away, and related to it by complex conjugation, while for every 
red element there is a map to the element with row number a distance +1 away 
and column number a distance —1 away, and related to it by complex conjugation 
times —1. Assuming that it is possible to use the same map a to increase the column 
number for each block of the matrix it follows, in addition to (4.1.7), that 




(4.1.22) 


for a blue matrix cell (/i, u) and 




(4.1.23) 
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for a red matrix cell (/r', u') where is the action of the map a on the m-body state 
l/r). Hence a blue cell corresponds to those given by (2.1.37) with 


= {Tfi\Hk\Tuy (4.1.24) 

for some /i, u while a red cell corresponding to (2.1.38) can be written as 

{y\Hk\Tu') = -{Ty\Hk\iy'y (4.1.25) 


for some p', z/'. Setting p' = Tp' and comparing equations (4.1.22), (4.1.23) with 
(4.1.24) and (4.1.25) gives 

T\p) = (4.1.26) 


Equations (4.1.22) and (4.1.23) imply that the random variables of the potential are 
subject to the condition 


^jkik 


(4.1.27) 

when the condition is yielded by applying (4.1.22) and are subject to the condition 




= -v„ 


"Ok 




(4.1.28) 


when the condition is yielded by applying (4.1.23). Hence for the eGSE the random 
variable Vj^i^ does not satisfy the same condition for all jk, ik- This distinctive feature 
of the eGSE indicates that it is necessary to count the number of sets of {jk,ik} 
for which the two conditions (4.1.27) and (4.1.28) overlap in order to calculate the 
moments of the level density of the eGSE. The overlapping of these conditions gives 
the symplectic zeros’’’’ which are the values of the random variable {jk-, ik] for which 
Vjki^ = implying = 0. In other words, it may be that the sets of k 

single-particle states j,i are present in m-body states defining a blue cell as well as 
other distinct m-body states dehning a red cell. This situation only arises for k <m. 
Eor k = m there can be no overlap between the conditions (4.1.27) and (4.1.28). 
When making the transition from studying the eGUE to studying the eGSE the hrst 
difference, and the most signihcant hurdle, is this overlap in the conditions on the 
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random variables 

Permutations 

As shown in section 2.1.2 the assumptions underlying the symplectic symmetry of 
H coincides with a fermionic time-reversal invariant quantum system. If there are 
m fermions in the system the state space will consist of m-body fermionic states. 
Dehning the system by a fc-body potential with I non-degenerate single-particle states 
there must be non-degenerate states in the basis describing the system and the 
matrix elements of the symplectic potential are given by := {fi\H\u). Since K is 
even (see section 2.1.2 ) it follows that H is even as well. If / = 2n is even and m 
odd, Glaisher’s Rule gives 

( l\ jo mod 2 if Z is even and m is odd; 

\m) y ([ 2 / 2 j) ^ otherwise. 

which ensures that there are an even number of m-body states and therefore an 
even number of rows and columns in the matrix representation of the potential H. 
This is just an example of the interesting dependency of m, Z and that occurs 
for potentials of the eGSE. It is known that a basis can be chosen in which the time 
reversal operator K takes the form of (2.1.28). With K taking this block diagonal 
form, the operator T = KC dehning time-reversal invariance, in addition to taking 
the complex conjugate of a state, maps one basis vector to another and multiplies 
it by -1. The operator being that which applies time reversal twice to a state, 
maps a basis vector to itself times a phase of -1 sans complex conjugation (see section 
2.1.2). By letting the state labels be the set {1,2,...,Z} for the basis in which H 
takes the canonical symplectic form (2.1.33) it must be the case that the time reversal 
operator T maps any given m-body state to another, times at most a phase. Hence 
the action T|/r) = is limited to permuting the m single-particle labels in the set fi 
and multiplying the result by some phase Assuming that this phase is the same 
for all states it follows that 


(4.1.29) 


T/J, = 


e*‘^cr(/i) 


(4.1.30) 
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where (j(-) is a permutation on the set {1,2,..., 1}. Applying = —1 gives 




SO that 


with p = 1, 2, 3, etc., and also 


2p- 1 


-TT 


a^ = l 


(4.1.31) 

(4.1.32) 

(4.1.33) 


which implies that the permutation map a is a pairwise permutation map. These 
permutations can be thonght of as single-particle state maps which exchange one 
single-particle state with another and vice versa as a function only of the initial state. 
For example, the pairwise permntation map 


a = 


1 2 3 4 ... 2n - 1 2n 

2 1 4 3 ... 2n 2n-l 


(4.1.34) 


exchanges the single-particle state 1 with 2, 2 with 1, 3 with 4, 4 with 3 and so on. 


The Zeros of the Random Variable 

The considerations so far for the second moments of the symplectic ensemble suggest 
that the eGSE involves far more snbtlety than the eGUE and eGOE. A snmmary 
of the reqnired approach to problems relying on the second moments, inclnding all 
calculations of the higher moments of the eGSE, can be expressed as follows: 


Fix T to canonical form 

-)■ 

Define a map a 

-)■ 

Gount Zeros 

-)■ 

Galculate Moments 


The form of T = KC is hxed in order to predetermine the form of the symmetry 
relations on H, namely (4.1.24) and (4.1.25) which are the analogue of the symmetry 
conditions of (2.1.33) written in terms of T. Next a map a is chosen to determine 
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the action of T. Any resulting zeros from an overlap of (4.1.27) and (4.1.28) need to 
be identified to completely characterise the second moments VjjVjE for all j, i,j', i'. 
Hence a “zeros theorem” needs to be stated for each permutation map a in order to 
identify the fc-sets {j,i} for which the random variable is not in fact random but 
instead equal to zero. Given these additional steps for the eGSE however, the form 
of the equations do not change dramatically. As for the eGUE, when Vji and are 
not related by (4.1.27) or (4.1.28) 

VjiVj'i' = ^ / VjiVj>->exp exp dvjidv->-' = 0 (4.1.35) 

since the gaussians are even functions, making the total argument of the integral odd 
with respect to the integration variables. Likewise if Vji = then Vji and n^-q' are 
complex and their product yields two even functions of opposite sign, again giving a 
zero average. However, for 

Vji = ±v*,. (4.1.36) 

the average becomes 

'Vji = j v%ex^ ^ (4.1.37) 

The sign of (4.1.36) is determined by the permutation map a. By substituting 
r|/i) = e^‘^a{ij) from (4.1.30) the socalled blue matrix elements with symmetry 
condition given by (4.1.24) obey 

itilHW) = (e‘*<Tq)|ff|e**<T(0>- = {a(y:)\H\a(v)Y (4.1.38) 

while red matrix elements (4.1.25) obey 

(4.1.39) 


implying 


{a{^^)\H\v) = -(/i|i7|a(z/))L 


(4.1.40) 
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From these it follows that 

(4.1.41) 

for fc-tuples j,i which are elements of the m-particle states determining blue matrix 
elements and 

"ji = (4-1-42) 

for fc-tuples j,i which are elements of the m-particle states determining red matrix 
elements. This finally yields the non-zero second moments as 

VjiVjTi^ = sgn(ji) ■ (4.1.43) 

Hence by the above argument and by (4.1.27) and (4.1.28) it can be concluded that 

Vjivy^ = 5-5- -t- 5^34 ■ sgn(ji) ■ 5^.^(y)5.^(./) (4.1.44) 

where the amplitude of 5^34 gives the purely Symplectic component of the average. 
The sgn function determines the sign of the /c-tuple, where the sign is positive for 
tuples resulting from (4.1.41), negative for those resulting from (4.1.42) and zero for 
those tuples satisfying both. 

4.2 The Embedded GUE 

4.2.1 Introduction 

Given the results from sections 4.1.2, 4.1.3 and 4.1.4 it follows that the second 
moment for the embedded gaussian unitary, orthogonal and symplectic ensembles 
can be written as the compound expression 

+ <^54 ■ sgn(ji) ■ (4-2.1) 

This is just the linear combination of equations (4.1.16), (4.1.21) and (4.1.44). The 
second moment is therefore greatly simplihed in the case of the eCUE, reducing to 
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read just 



(4.2.2) 


It is for this reason that calculating the moments of the embedded gaussian unitary 
ensemble is the least difficult of the three cases. This thesis will therefore focus 
primarily on the eGUE case, even though the methods introduced can be applied to 
other symmetry classes as well. To begin with, a known result for the normalised 
fourth moment, or kurtosis, of the level density for the eGUE will be derived using a 
simpler alternate calculation to that found in [BRWOla]. This will serve as a platform 
for introducing both the basic definitions and the new methodology involving particle 
diagrams which this thesis aims to explain. In subsequent sections the same methods 
will be used to calculate the sixth and eighth moments as well, albeit with the addition 
of some extra methodological features. Eor the case of the eGUE we will consider 
m spinless fermions in a system with m -C I single-particle states, all interacting 
under the action of a A:-body potential {k < m) with an identical gaussian p.d.f 
determining its independent entries. The single-particle creation and annihilation 
operators are a] and aj respectively with j = 1,... ,1. As introduced in prior chapters, 
a useful shorthand notation for products of these is the abbreviation j = (ji,..., j^), 
aj = ajf, ... aj^ (similarly for i). A corollary of this is 



(4.2.3) 


For the sake of continuity and comparison with extant literature, notational traditions 
will be preserved for the calculation of the kurtosis by writing the orthonormal m- 


particle states as \fi), |z/), |p), etc. where each state takes the form ... u^jO) with 


|0) denoting the vacuum state and the restriction 1 < ji < J 2 < • • • < im < k The 
fc-body potential is given by 




■ ■ ■ ^jk^^k ■ ■ 


. a 7 a,- 




(4.2.4) 
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which will be abbreviated to 

= (4.2.5) 

as before. In this case all the m-particle states are assumed to contain only non¬ 
repeating single-particle states and since the number of allowed single-particle states 
is /, the dimension of the m-body state space is Y = . For the embedded CUE 

ensemble presently under consideration the only symmetry condition on the potential 
is that it be hermitian {fi\Hk\i') = {u\Hk\fi)* for all fi, v giving the simplihed second 
moment = Vij derived in section 4.1.2. As in canonical RMT (hxed k = m) the 
matrix elements not related by hermitian symmetry are uncorrelated i.i.d complex 
gaussian random variables with mean zero and variance = 1 (wlog). 

4.2.2 Particle Diagrams 

Of particular use in calculations is the abbreviation 

■= {fi\Hk\a){p\Hk\u) = {p\a]ai\a){p\alaj\u). (4.2.6) 

Here summation over the repeated indices i,j is implied. For to be non¬ 

vanishing |(j) and Ip) must both contain the k single-particle states included in i, and 
Ip) and |z/) must both contain the k single-particle states included in j. In addition 
ai\p) and aj\a) have to contain the same single-particle states implying that |p) and 
|cr) coincide in the m — k single-particle states not included in i or j. The same 
applies to |p) and \u). These relations are illustrated in Fig. 4.1a where solid bonds 

- connect states sharing m — k single-particle states and dashed bonds- 

connect many-particle states sharing k single-particle states. Note that in this figure 
the overlaps indicated by neighboring bonds are disjoint, e.g., the single-particle states 

of i form the overlap a - -p but are excluded from the overlap p- a because 

the state |(j) contains m non-repeated single-particle states. The “particle diagrams” 
drawn in this way form an essential ingredient for evaluating the moments of the 
level density. With the odd moments being zero the normalised 2n’th moments of 
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the level density are given by 


SO the fonrth moment, or knrtosis, which will be calculated next is 



(4.2.7) 


Wtrm 


(4.2.8) 


The terms in the nu m erator and in the denominator of k can be evaluated using 
Wick’s theorem which expresses the Gaussian average as a sum over all ways to draw 
contraction lines between the factors 14. This gives 


" 2 ', _ 


n 

(trV.Vk) 


tr(V^) = (trl4l4l4l4) + (trl4f4l4l4) + (trl4l4l4l4). 


_ n n 


n 


(4.2.9) 

(4.2.10) 


Here the brackets (.) are used as an alternative notation for the ensemble average. 
Note that the hrst two contributions to (4.2.10) coincide due to cyclic invariance 
of the trace. Wick’s theorem then allows us to compute the averages of contracted 
matrix elements as if the remaining elements were absent. The result can be expressed 
in terms of (4.2.6), leading to 

tr(7l!) = ^ ( 4 . 2 . 11 ) 

in the denominator, while in the numerator 

tr(77^) T ^^-pypuAupyp (4.2.12) 


with the summations over repeated indices /i, p, a implicit. One can alternatively 
derive (4.2.12) from hrst principles by observing that the random variables Vji are 
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gaussian so ^ and / x^e~°‘^^dx = Given the restrictions 

from Fig. 4.1a it can be seen that it is necessary to sum over all |/i) and \p) sharing 
m — k single-particle states in order to calculate (4.2.11). There are N = states 
in the sum over all possible |/i), ways to choose the overlap with \p), and 

ways to choose the rest of \p). Hence the result is 


The trace in the numerator of k is given by (4.2.12). The calculation for the hrst 
term is almost identical to the calculation yielding (4.2.13) giving 


9 4 4 

^^aapp-^aafifi 



^ — m + 


(4.2.14) 


(summation implicit) so that the only remaining term needed to complete the calcu¬ 
lation for K is This term, as well as subsequent quotients dehning the 

sixth and eighth moments, requires the summation of a series of binomial expressions 
not all of which are simple enough to write down, as has been done with (4.2.13) and 
(4.2.14). 


4.2.3 Binomial Arguments 

To determine which of these binomial expressions survive in the limit of large I a 
new quantity called the argument will now be dehned as the power of a binomial 
expression in I in the limit / —)■ cx). For a quotient to give a non-vanishing result in 
that limit the argument of the numerator must be at least as large as the argument of 
the denominator. The argument can be obtained using Stirling’s formula. Taking the 
dimension of the state space N = as an example and applying Stirling’s formula 
limn^oonl = y/27m(^^^ the value of the argument of N is taken as the power of 
N{m\) ~ is so arg(Y) = m. More generally 
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Figure 4.1: (a) illustrates the particle diagram implied by = (/i|ajaj|(j)(p|ataj|i/) 

of (4.2.6) and (4.2.12), and (b) illustrates the particle diagram implied by the factor 
in (4.2.12). Each bond between compound m-body states represents a set of 
single-particles shared by both states. 


arg 





n 


or alternatively in terms of multinomials 


arg 


/ I — a 
Ml ^2 h ■ 




(4.2.15) 


(4.2.16) 


Note that any additional factors independent of I do not have any impact on the 
argument - they contribute only hnite coefficients to the product and therefore do 
not affect the power of /. Erom (4.2.13) it then follows that arg[(^tr(V)?))^] = 2k. 
Since it is expected that the value of the fourth moment will converge to some hnite 
value the term ^tr(V)f) should have terms with arguments equal to 2k and possibly 
some terms with an argument less than 2k, which indeed turns out to be the case. 
Those terms with an argument less than 2k will be ignored, contributing values of 
order no higher than to k as / —)■ cx), whereas those with an argument equal to 2k 
must be calculated. Using (4.2.14) one easily sees that after division by N the term 
has the argument 2k as expected. 

Eor the second summand Api,p„A„pi,p in (4.2.12) non-zero contributions arise if 
the indices of A^^p^ obey the same restrictions as introduced earlier with reference 
to (4.2.6) and depicted in Eig. 4.1a. The analogous restrictions for non-zero A^pup 
are given by Eig. 4.1b. Together Eigs. 4.1a and 4.1b form the particle diagram for 
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Afj,,ypaAaij,up, formally defined as the collection of bonds incorporating all restrictions 
for the indices. Since there will be occasion to use the combined diagram in later 
calculations for the higher order moments as well, it will be referred to subsequently 
as the standard diagram. Of interest is the number of unique m-body states satisfying 
this diagram (giving the value of the trace) but only the term for which the argument 
of this number reaches its maximal value 2k is needed. This is the term for which 
the maximal number of single-particle states can be chosen which satisfy all of the 
bonds of the particle diagram, ie. the term for which there are maximal degrees of 
freedom in the choice of single-particle states which satisfy the bonds of the diagram. 
From Fig. 4.1a it can be observed that the state |/r) determines the k single-particle 

states in the state \i') which are not determined by the bond v - p so that \p) and 

Ip) together fully determine the single-particle states contained in \i>), and similarly 
for the state ja). Hence all the single-particle states satisfying an instance of the 
diagram are completely determined by \p) and \p) and the argument of the diagram is 
maximal when \p) determines the smallest possible number of single-particle states in 

Ip) and vice versa. This means that e.g. the bond p- u of Fig. 4.1a must share the 

minimal number of states possible with the bond u - p of Fig. 4.1b, and likewise the 

bond p-cr of Fig. 4.1a must share the minimal number of states possible with the 

bond cr- p of Fig. 4.1b. The minimum overlap possible between the smaller pair is 

zero, while the minimum overlap between the larger is s = |/c — (m — k)\ = m — 2r 
where r := min(/c,m — k). This overlap being minimal additionally implies that the 
larger of the two bonds between two nodes (states) contains all single-particle states 
included in the smaller one, leaving r states participating in both bonds and s states 
participating only in the larger of the two bonds. Analogous reasoning applies to each 
of the corresponding pairs of bonds of Fig. 4.1a and Fig. 4.1b. The sets of s states for 
each pair must be disjoint relative to the sets of r overlapping states, as a consequence 
of neighboring bonds in the same diagram 4.1a or 4.1b being disjoint. Hence the 
possible choices for states are given by partitioning the I available single-particle 
states into one set of s states shared by both |p) and |p), and four sets of r states. 
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Thus to leading order the expression can be written as the multinomial 


A A rsj 


I 


\ S T' ^ T' } 

l \ fl — m\ fl — m — r 
ml \ r j\ r 



(4.2.17) 


Recalling N = the argument of the corresponding term in the numerator of 


K 


is then arg ^ Y, = 2r = 2 ■ min{k, m — k) so it is only ior k < m — k 

that the argument of this term is equal to 2fc, while for k > m — k it is always less. 
Eitting all the surviving terms into the expression for k gives the limit form of the 
fourth moment as / —)■ cxo 


lim K, 

N—^OO 


T 4 A 

TV—^CXD 1 A 

V 4^ 





(4.2.18) 


which corroborates the result found by Benet et. al. [BRWOla] using the eigenvector 
expansion method and for 2k > m agrees with what is expected using the method of 
supersymmetry, namely k = 2. 
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Figure 4.2: k against k/m for m = 6,8,10,12 showing how the kurtosis has a semi¬ 
circular value, K = 2 for k/m > \ after which it transitions to a gaussian kurtosis, 
K = 3 at fc = 0. Higher values of m give faster convergence to the semi-circular 
moment. 

4.2.4 Paths and Loops 

Particle diagrams represent diagramatically the factors of (4.2.6) which form 

the constituent parts of (see (4.2.12) for instance). What will be needed next 
in order to calculate higher moments for embedded ensembles is an extension of 
the notation of particle diagrams to combine arguments with particle diagrams in a 
mathematically robust way. In so doing, it becomes possible to calculate any diagram, 
no matter how complex. This will be done by using arguments (powers of 1) and 
diagrams (representations of as before, with the additional concept of closed 
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loops. It is a necessary to introduce closed loops on diagrams because as we look 
to calculate higher and higher moments the resulting diagrams become increasingly 
complex and less responsive to such a simple analysis as that which yielded the fourth 
moment. Closed loops form the unbreakable link between a diagram and it’s maximal 
argument. 

As explained in section 4.2.2 a particle diagram is a set of bonds drawn between 
m-body states indicating the minimal overlap between them. Each m-body state in 
the diagram can be thought of as a node. In calculations of the moments of the level 
density each node is connected by four bonds, two of size k and two of size m — k. 
See Fig 4.3 for a representation. This is a property that emerges from the application 
of Wick’s theorem to which can always be expressed in terms of sums of factors 

of = {p.\H\a){p\H\h'). Since from Erst principles the trace is simply 

tri 72 n = Ha^^a^Haia^Ha^ai ■ ■ ■ Ha2„ai (4.2.19) 

every index appears twice, which directly implies that each index will appear twice 
in the resulting expansion of (4.2.19) in terms of Hence apart from a single 

exception (discussed next) every index will appear once in one factor of A and a 
second time in another distinct factor of A. This in turn implies that every node in 
every diagram will take the form of Fig 4.3 wherein the node (a set of m state labels) 

is completely determined by a fc set-disjoint from the set- oi m — k labels. 

That these must be disjoint is given by the fact that the nodes represent fermionic 
states which are sets of m non-repeated labels. Morever, since there are two such 
disjoint sets connected to each node, one for each appearance of a state as an index 
of A, these must overlap entirely in order to conserve the number of states in the 
node as m. In other words, with reference to Fig 4.3, the state \fi) is the union of the 
two disjoint black bonds 

^ = -u- (4.2.20) 

from the hrst factor of A containing p as an index, as well as the union of the two 
disjoint blue bonds 

p = - 


— U 


(4.2.21) 
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from the second factor of A containing the state fi as an index. Hence, for example, 

the m — k states determined by the blue solid bond - must themselves form a 

disjoint set; those which are shared with the black solid bond- fi and those shared 

with the black dashed bond-In this way, the single-particle states which form 

the set denoted by the blue solid bond will “travel” (read: be a member of) through 
either the black dashed bond or the black solid bond, and onwards to the next node 
(see Fig 4.4). At the next node this process will be repeated and the single-particle 
states will seperate again into disjoint sets and travel further, continuing along a path 
around the diagram until forming a loop. 

For the singular case where a single factor of A^^pp includes both instances of a 

label the two bonds-connecting fi and p are equal and are represented by a so 

called “tail” p - p, and the fc-bonds-are implicit in that they connect each 

state with itself, as opposed to the more common case where for example p shares a 
fc-bond with two distinct states, neither of which is equal to p. 

These loops on the particle diagram indicate the single-particle states shared by 
all bonds and nodes (states) which participate in the loop. There are only ever a 
hnite number of possible loops in any diagram. This relates directly to the argument 
of the diagram, since the argument counts the degrees of freedom in a diagam, i.e. the 
number of labels which together determine, through bonds, all states in the diagram. 



Figure 4.3; Each node in a particle diagram takes the basic form shown. Since the 
blue bonds are disjoint and together completely determine the m single-particle labels 
in \p) they must overlap completely with the set dehned by the black bonds, which 
are also disjoint and completely determine the single-particle states in \p). 
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Hence the number of elements in a loop is the number which that loop contributes to 
the total argument. In this way loops will form a fundamental tool for (a) expressing 
the argument as a sum of loops and (b) maximising the arguments by optimising the 
number of elements in each loop. 

4.2.5 Example: The Fourth Moment with Loops 

To illustrate the process of calculating arbitrarily high moments using arguments, 
particle diagrams and loops, an example will now be made with of fourth moment, 
which has already been calculated from first principles as (4.2.18). It has been 
seen that the fourth moment, or kurtosis, k, = ^tr(iif|)/(|^tr(Tr|)^ where in the 
denominator one has 

= ( 4 - 2 . 22 ) 

and in the numerator 


— ‘^^crapp^acriiii 


+ 4 


A 

fiupa^a^up' 


(4.2.23) 



Figure 4.4: The fact that each node in a diagram takes the form of Fig 4.3 implies that 
the combined diagram can always be “flattened” into the form of a cyclic lattice. This 
implies that every single-particle state label in the particle diagram is the member 
of a loop which travels around the lattice along those paths not excluded by the 
fermionic nature of the nodes (states). 
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Figure 4.5: The particle diagram of the term is the standard diagram, 

as seen before in Fig 4.1a and 4.1b but now with both factors illustrated as a single 
combined particle diagram. 


As forecast in the previous section (i) both components of the fourth moment can 
be expressed entirely in terms of the quantity A^^p^^ and moreover (ii) every index 
appears twice, either as an instance of a tail with repeated indices of the form Aaabb 
or once each in two distinct factors of A such as the second term Apyp„A„pyp of 
(4.2.23). For terms of the form Aaabb no further sophistication is necessary because 
these simply add a single factor of cases where only 6 is a free choice, 

or alternatively cases such as (4.2.13) where both a and b have not 

yet been determined by other bonds in the diagram. These considerations give 



(4.2.24) 


and 



(4.2.25) 


The term Ap^p^A^pup can be expressed as a particle diagram as in Fig 4.1a and 4.1b, 
however this time the diagram will be drawn as a single combined picture. Fig 4.5. 
This is the form in which it is possible to calculate the moments by identifying all 
loops within the diagram. Because of the fermionic nature of the states, loops cannot 
contain two bonds from the same factor of A, but apart from this caveat all loops are 
permitted. The loops for Fig 4.5 are illustrated in Fig 4.6. Lines curving outwards 
indicate a path through a bond of the red component A^pup while curves bending 
inwards indicate that the path is crossing a bond in the blue component Ap^p^. The 
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(a) The loop JIdfi defin¬ 
ing the set of single-particle 
states shared by the bonds 

yU _ a and jj, _ a of 

Fig 4.5. 


(b) The loop ap^ defin¬ 
ing the set of single-particle 
states shared by the bonds 

a _ p and a _ p of 

Fig 4.5. 


(c) The loop Jlufi, defin¬ 
ing the set of single-particle 
states shared by the bonds 

p _ u and p _ u of 

Fig 4.5. 



(d) The loop up^ defin¬ 
ing the set of single-particle 
states shared by the bonds 

u _ p and V _ p of 

Fig 4.5. 



(e) The first loop papufi 
defining the set of single¬ 
particle states shared by the 

bonds p _ cr, ex _ p, 

p _ v and z/ _ p of 

Fig 4.5. 



(f) The second loop papufi 
defining the set of single¬ 
particle states shared by the 

bonds p _ (J, a _ p, 

p _ u and V _ p of 

Fig 4.5. 


Figure 4.6: Illustration of the loops of the factor Afj^i,puAufj,up seen in (4.2.12) with 
the factor A^^per being represented by the particle diagram of Fig 4.1a and the factor 
Aapup being represented by the particle diagram of Fig 4.1b. Each bond between 
compound states represents a set of single-particle states shared by both of the 
compound states, and each loop represents the set of single-particle states shared by 
all of the compound quantum states which form the nodes of the loop. 
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sum of the single-particle labels in every loop in the diagram is necessarily the sum 
of all dinstinct single-particle labels contained in the diagram and hence is equal to 
the argument. Labelling each loop gives 

a = [loIl 
(5 = ap^ 

7 = Jwfi 

6 = vpl> 

e = ]Io~p0 

K = papp^ (4.2.26) 

so that 

arg[A^^p(jA(j^^p) = a-\-l5-\-'y-\-5-\-e-\-K. (4.2.27) 

Since the total number of single-particle states in each bond-is k and the sum 

of all single-particle labels in each of the solid bonds-must be m — /c it becomes 

possible to write the following “conservation equations” wherein the sum of the 
elements (read single-particle labels) in all loops passing through a given bond must 
equal the total number of elements in the bond (either k or m — k). Now each loop 
is taken to denote the number of elements in that loop. The conservation equations 
for the bonds of the red component of the diagram read 


a + e = m — k 

(4.2.28) 

13 + n = k 

(4.2.29) 

-y + K, = k 

(4.2.30) 

5 + e = m — k 

(4.2.31) 


while the conservation equations for bonds belonging to the blue component of the 
diagram are 


a + K = k 


(4.2.32) 
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P^e = m-k (4.2.33) 

7 + e = m — fc (4.2.34) 

5 + K = k. (4.2.35) 

Plugging equations (4.2.28) and (4.2.29) into (4.2.27) gives 

arg{A^^p^A^^^p) = m + 7 + 5. (4.2.36) 

Now additionally substituting (4.2.30) and (4.2.35) gives 

arg{Ap„p„A„p„p) = m + 2k - 2 k (4.2.37) 


which is maximised when k = 0 , implying that to attain the maximum argument 


k — K = k 

(4.2.38) 

k — K = k 

(4.2.39) 

k — K = k 

(4.2.40) 

k — K = k 

(4.2.41) 

m — k — 'y = m — 2k. 

(4.2.42) 


The hnal step is simply selecting the single-particle states for each loop from the total 
number of possible single-particle energy levels I to give the leading order term 


A A 

■^fj^upa-^a^up 


^ I ^ 
,m — 2k k k k k j 


(4.2.43) 


This shows, as found already by direct arguments, that the fourth moment is 


lim K 

N^OO 


oT A A 

^ ^ -^aapp-^acTpip 


+ 


A 


2 — 


( A A 


+1( 


I 

m—2k k k k k 


(t(.9(T)pr‘))' 
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= 2 + 



(4.2.44) 


The above explanation illustrates the complete application of the method of particle 
diagrams. This involves firstly the notion of an argument which identihes the order of 
magnitude of terms in powers of /. Numerators in a quotient must necessarily attain 
the same argument as the denominator or they will not survive in the limit / —)• oo. 
Those terms which will not survive can be ignored. Secondly, the method of particle 
diagrams maps every term into the form of a diagram of nodes and bonds. 

Finally, the method of particle diagrams is completed by identifying all loops 
within a diagram and maximising the argument by maximising the sum of the (number 
of single-particle labels in) loops as seen in the example of (4.2.27). After this the 
practitioner must simply select the components of each loop from the set of possible 
single-particle state labels {1,2,..., /} as seen in the example of (4.2.43). 





Chapter 5 

The 6’th and 8’th Moments 


5.1 Sixth Moment 


If the fourth moment conhrmed that the method of particle diagrams can greatly 
simplify calculations, the sixth moment illustrates the flexibility of the method as well 
as its ability to scale. For the sixth moment a transition is made from flat particle 
diagrams to three-dimensional graph-like diagrams. Indeed, the terminology used 
will be increasingly that of graphs and sets. The m-body states will be nodes on 
the graph. Graphs will contain paths - seqnences of neighbouring nodes. And the 
concept of loops on the particle diagrams, where the first state in a path equals the 
last state, will become an increasingly important tool for maximising the argument 
of attendant binomial expressions as the complexity of the particle diagrams become 
more complex. The sixth moment of the level density is given by 


Using Wick’s theorem as before furnishes 


(5.1.1) 


n n n 


n 


n 


n 


tr(U,6) = 2(trWWWW14W) + 3(trW14WW14W) + 6(trWWWWWW) 
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+ 3{tTVkVkVkVkVkVk) + {tiVkVkVkVkVkVk). (5.1.2) 


Here the prefactors indicate the number of equivalent diagrams that can be obtained 
by cyclic permutation of the trace. Written in terms of the summands are 




9 4 4 4 

^ ^ptqq^tvuu-^vpww 


+ 3A 


ptqq-^uwvv 


i^-^tpwu) 


{.-^twvu-^upwv^ “ 1 “ ^-^putq-^qwvt-^upwv 
“t“ ^pvuq-^qwvt^tpwu 


(5.1.3) 


Terms involving A’s with identical first and second (or third and fourth) indices 
simplify as they give a contribution only if the two other indices coincide as well. 
For instance for Apt^g = {p\(i](ii\q){qlalajlt) to be nonzero the states aj\p) and aj\t) 
both have to coincide with ai\q). Adding the single-particle states with indices in j 
then gives coinciding \p) and |t). Using this idea as well as the reasoning leading to 
(4.2.13) the first two terms in (5.1.3) can be evaluated to give 

2 ApiqqA^yiiuAyp^^ -|- 3 Ap^qqAuwvv 


i^Afp^y^ 5 


m, 


'm\ fl — m + 

ifcjl k 


(5.1.4) 


The term SAptqq (AtwvuAupwv) 

For the third component QAptqq (AtwyuAupwv) it is similarly required that \p) = \t) 
for a non-zero contribution. The particle diagram for this term is illustrated in Fig. 
5.1a. Note that it is nearly identical to the particle diagram of Fig 4.5 except for the 

addition of a tail p - q which adds a factor obtained as for Eq. (4.2.13). 

This additional factor added to the expression we know already gives the leading 
order expression when / —)■ cx) as 


Ap,,, ~(j.j( f. j(jrrrr, 


(5.1.5) 




















5.1. SIXTH MOMENT 


93 




Figure 5.1: (a) illustrates the particle diagram for the term QAptqq {Aty^^uAupwv) ■ This 
is almost identical to Fig. 4.1 but with the equivalent of Fig. 4.1a and Fig. 4.1b 
jnxtaposed on the same diagram as in Fig 4.5. The additional tail adds a factor 


0 ( 


k 


(b) illnstrates the particle diagram of SAputqAqwvtA 


upwv 


where all the 


bonds implied by the expression are illustrated in a single 3D triangular prism. Using 
the diagram it becomes simpler to identify the single-particle states which mnst 
overlap maximally in order to maximise the argument of the whole sum. 


The argument of this term is A: -|- 2r where r = min(k,m — k) is dehned as before, so 
it will only survive in the limit of h as / —)■ oo for k < m — k. 

The term AputqAq^vt.^ upwv 

The particle diagram for the fonrth term AputqAqu,vtAupwv of (5.1.3) is illnstrated in 
Fig. 5.1b as a single three dimensional triangular prism denoting the interrelated 
conditions on the states that mnst be satished for the expression to be non-zero. 
The three faces correspond to the three factors. Since these are fermionic states as 
before, adjacent bonds on the same face cannot share single-particle states. In order 
to maximise the number of participating single-particle states in the sum, and hence 
the argument, the four overlaps of k states between |u), \t), \q) and |m) are chosen to 
be disjoint. This imperative will be shown later using the full method involving loops, 
but for now it will be momentarily assumed in order to complete the calculation of the 
result sans loops. Next choose the m —2k additional states participating in the bond 

V - w bnt not in v - w. Note that the many-particle state |u) is given by the 

combination of the bonds v - w and v -and equivalently by the combination 
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of V - w and v - u. These combinations can coincide only if v - u involves the 

overlap v - 1 as well as the aforementioned m — 2k states. Hence the m — 2k states 

participate not only in |r;) and |tc) but also in |m). An analogous argument shows 
that they also have to be included in \p), and thus in all many-particle states of the 
“left” face in Eig. 5.1b. To fully determine |m) and \p) it is still necessary to choose 

the k states participating in the “left” of the two bonds u - -p, which altogether 

gives m — k = m — 2k+ k states in addition to the original four sets of k states. 
Considering the “right” face in an analogous way one obtains the same choice of 
m — k states but now broken down differently into a choice oi m — 2k and a choice of 
k states. Hence our choice is restricted to selecting four sets of k states and one set 
oim — k states from the I available states, and then splitting the latter into sets of 
m — 2k and k states in two independent ways. This gives 

How to calculate this utilising loops? As before the hrst step is simply to identify all 
possible loops in the diagram, in this case the prism-like diagram of Fig 5.1b. These 
loops are illustrated in Fig 5.2. Giving each loop a label indicating the number of 
single-particle state labels contained in the loop gives 


a = utqpu 
(3 = uvwp^ 

7 = vutqpwv 
6 = vutv 
e = qpw^ 

K = vtqwv 


p = qputvwq 
0 = wpuvtqw 
p = ^t 
71 = vwh 
u = pu^ 

(5.1.7) 


Since the number of labels passing through any bond-must sum to k and the 

number of labels passing through any solid bond-must sum to m — fc it is 

straightforward to write down the following conservation equations which restrict the 
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(d) The loop vutv (e) The loop qpw^ (f) The loop vtqwv 




(j) The loop vwi (k) The loop pu^ 

Figure 5.2: Illustration of all possible loops within the prism-like particle diagram of 
Fig 5.1b corresponding to the term AputqAq^yt^upwv 
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values of the loops to 

a + rj + 00 = k (5.1.8) 

a + 'j + 6 + ri = m —k (5.1.9) 

a + 7 + (p = m —fc (5.1.10) 

a + 'y + e + r] = m —k (5.1.11) 

/3 + (j) + 00 = k (5.1.12) 

(3 + 'y + 5 + (j) = m —k (5.1.13) 

(3 + 'j + n = m —k (5.1.14) 

/3 + 7 + e + 0 = m —fc (5.1.15) 

e + K + ri + (j) = k. (5.1.16) 

The argument of the diagram is simply the sum of all degrees of freedom, which is 
the sum of the number of labels in each loop 

0‘^9^-^putq-^qwvt-^upwv') Q? + /5 + 7 + (5 + e + K + ?7 + 0 + (p + 7r + u;. (5.1.17) 

Substituting (5.1.9) and (5.1.12) yields 

9 (^-^putq-^qwvt-^upwv^ T^T^cT^pT 71 (5.1.18) 

and hnally inserting (5.1.16) gives 

0>7'9^-^putq-^qwvt-^upwv) k ^ 77 Tj (j) (5.1.19) 

SO that the maximal argument m + 3/c is achieved when ip = tt = k and 7 = 0 = 0. 
Plugging these values back into (5.1.32 - 5.1.16) reveals that ip = 7T = 5 = e = k and 
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K = rj = (j) = 0. Additionally 


a + 'j = m — 2k 
+ 7 = m — 2/c 
a + u = k 

^ + u = k (5.1.20) 


so that a, jS and 7 are functions of u and 


2^putq2^qwvt2^upwv ~ ^ ) 


= E 


= E 


^k k k k oj a{u) (3{u) 7 ( 0 ;) / 

' ) 

Jikkkuk — uk — um — 3k+ uj 

I \/m — k\ /k\ km — 2k'' 

^kkkkm — kj\ k J \oj J \ k — u ^ 

2 


^ / \ m-k 

,k k k k m — k) \ k 


(5.1.21) 


where the hnal step uses Vandermonde’s identity = (t) ^ (E 

the same result claimed earlier in (5.1.6). 


The term ApvuqAq^wtAtpwu 

For the final term ApyuqAqwvtAtpwu of (5.1.3) illustrated in Fig. 5.3 it is useful first to 
notice that the argument never exceeds 2m. This can be observed on sight since the 
states |tc) and |f) together define all the single-particle states in |m), |n), \p) and \q). 
And with an argument never exceeding 2m it can be concluded that this term will only 
contribute to the limit value of h for 3k < m. Since Iw) and \t) together determine all 
the single-particle states in the diagram and hence all the m-body compound states, 
there will be maximal degrees of freedom in the choice of the single-particle states 
in Iw) and |f) when |t(;) and |f) overlap minimally. This also determines when the 
argument for the entire sum (diagram) is maximal. Single-particle states shared by 
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both \w) and |t) must necessarily be manifested in the bonds which make up the 
diagram, so for instance if a single-particle state a is in both Iw) and |t), then it must 
have “travelled upwards” from the bottom node of the diagram at Iw), proceeded 

“through” (read: be a member of) the bonds w -p, w -f, w - q or w - u and 

continued towards the uppermost node of the diagram |f) “through” bonds q - 1, 

u - 1, p - -t or V - -t. More concretely, taking the bond w - p to represent the 

set oim — k single-particle states in both Iw) and \p) but not in |m), the bond w - u 

to represent the set of k single-particle states in both |tc) and \u) but not in \p), and 
so on, it can be seen that for any single-particle state a in both |tc) and |t) it follows 
both that 

a E w - pUw - V U w - qUw - u (5.1.22) 

and 

aEq - tUu - tUp - tUv - 1 (5.1.23) 

whereas for a given single-particle state d in \w) but not in |t) it follows both that 

a E w - pUw - V U w - qUw - u (5.1.24) 


and 

a^q - tUu - tUp - -tUv - -t. (5.1.25) 

This implies that states such as a in both \w) and |t) can be found in a connected 
non-repeating sequence of neighbouring nodes in the diagram starting at \w) passing 
through |t) and looping back to \w). Likewise for a single-particle state a which is 
in \w) but not in |t) one can trace a non self-crossing loop of bonds in the diagram 
starting from |t(;) and returning back to |tc) without passing through the node defined 
by the state |f). In conclusion, since it is imperative to maximise the number of 
single-particle states that are in |w) but not |f) and vice versa, and since the number 
of single-particle states which are an element of a loop excluding the node |t) is the 
complement of the number of loops passing through the states (nodes) Ire) and |f), of 
all the particle diagrams taking the form of Fig. 5.3 those with the largest argument 
are the ones consisting of the maximum possible number of loops from Ire) back to 
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|w) without passing through |t), and from |t) back to |t) without passing through 
|tc). Any diagram with even one less single-particle state than the maximal nnmber 
partaking in these loops will have an argument less than the diagrams containing 
the maximal number of loops, so these do not have to be calculated. As before 
it is prohibited to include a path containing bonds which are mutually exclusive; 
neighbouring bonds from the same factor of A. In Fig. 5.3 for example the bonds 

w - q and w - v are mutually exclusive as they both appear in the diagram due 

to the same factor in ApuuqAqwvtAtpwu- The k single-particle states in w - q 

are therefore the complement of the m — k single-particle states in w - v, the 

union being the set of single-particle states contained in |w). On the other hand 

the bonds w - q and w - p appear in the diagram dne to separate factors Ag^^t 

and Atpwu respectively so these are allowed to share single-particle states. Hence the 
loop wquvib is not a valid loop for the diagram of Fig. 5.3 as it contains mutually 



Figure 5.3: The particle diagram for the term Ap^uqAgwvtAtpwu which takes the form 
of a regular octahedron, or two square pyramids with a shared base determined by 
the plane on which the sub-diagram for Ap^uq is illustrated. The states \t) and jw) 
together determing the states |n) and \q) throngh the bonds dehned by Ag^yt just as 
they determine the states |m) and \p) through the bonds dehned by Atp^y Hence the 
maximal degrees of freedom in our choice of single-particle states from the set of I 
possibilities is bounded by 2m. 
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exclusive bonds, whereas wqpth is a valid loop. To summarise, a single-particle state 
which is an element of both |tc) and |t) must necessarily be an element of every bond 
in a loop passing through |tc) and \t). A single-particle state in |tc) and not in \t) 
must necessarily be an element of every bond in a loop from \w) to |w) not passing 
through \t) and conversely for a single-particle state in |t) but not in Since the 
argument of the diagram is maximised by maximising the number of single-particle 
states contained in every bond of a loop passing through only one of \w) or |t), the 
only diagrams which will contribute in the asymptotic regime are those where k 
single-particle states from each of the bonds in the loop 

Pi = wqp'ii (5.1.26) 


are equal to the k single-partical states in the bond w -g, each of the bonds in the 

loop 

P 2 = wuvih (5.1.27) 

contain k single-particle states identical to the single-particle states in the bond 

w - u, and so on for all the loops ps = wpvih with the bond p -u, p^ = tpqt 

with the bond t --p, ps = tuvt with the bond v - -t and pe = tuqt with the bond 

u - q. Borrowing set notation again one has for the case of (5.1.26) for example, 

that 

w - gflp - q = w - q (5.1.28) 


and 

w - qDw - p = w - q. (5.1.29) 

In this way the 3k single-particle states from |t(;) which participate in the loops pi, 
P 2 and p 3 do not overlap with \t) and the 3k single-particle states from \t) which 
participate in the loops p 4 , ps and pe do not overlap with \w). The remaining 
m — 3k single-particle states from Iw) and \t) which cannot be included in loops 
must necessarily be included in a loop through |tc) and \t). That is, the number of 
overlapping single-particle states between the m-body states |w) and |t) is m — 3k. 
This has profound implications for the value of the sixth moment and we will see a 
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related feature appear in the eighth moment as well. To make the comparison with 
the corresponding calculation (4.2.17) for the fourth moment more explicit, dehne 
r = min{k, m — fc) as before and s = m — 3r. We then have 


■^pvuq-^qwvt-^tpwu ~ 


^ ry' ry' ry ry 


(5.1.30) 


It should be noted that this time a relatively complex particle diagram has been 
evaluated with loops, albeit not by identifying all loops and maximising the argument 
as before, but instead by using arguments unique to the particular diagram at hand 
to determine what the size of each loop must be. It is useful to keep in mind that this 
is sometimes possible for larger particle diagrams, as it offers a shortcut to a result 
which would otherwise require a potentially far more laborious calculation. From 

(5.1.30) it follows that 


lim 

N—^oo 


T 4 A A 

]\j ^pvuq^qwvt^tpwu 




( m—k\ {m—2k\ 

k )\ k ) 



(5.1.31) 


Taking the quotient for h using the above expressions (5.1.4), (5.1.5) (5.1.21) and 

(5.1.31) gives the hnal result for the sixth moment 


lim h 

N—^OO 



(5.1.32) 
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Figure 5.4: h against k/m for m = 6,8, 10,12. Once again we see a domain k/m > \ 
for which the sixth moment takes the semi-circular value h = 5 and a transition 
thereafter towards a ganssian moment, h = 15 at k = 0. Higher valnes of m give 
faster convergence to the semi-circular moment. 

5.2 The Eighth Moment 

The eighth moment will exhibit the same features as the lower moments; a transition 
starting immediately after 2k = m from a semi-circular moment r = 14, to a gaussian 
moment r = 105 and a single term of tr taking the form 
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analogous to (4.2.17) and (5.1.30) and with s := m — = m — 4k. Of course we must 

again calculate products of ensemble averaged pairs of the matrix elements of 
for which there are now = 105 components in the sum. As may be expected, 

the method of particle diagrams can be applied analogously to the lower moments 
though with the diagrams now increasing in size and complexity. Some of these will 
be identical to particle diagrams which have already been seen in calculations of the 
lower moments, but with the addition of tails leading to extra factors of • 

This feature was seen in the transition from hgure 4.1 to hgure 5.1a for example. 
Eor progressively higher powers of the particle diagrams become progressively 
larger and, as for the diagrams illustrated in Fig 5.1b and Fig. 5.3 one can think of 
some of these as graphs in three dimensions, rather than as a flat structure. The 
larger structure of the subsequent diagrams makes it more difficult to identify just 
two states which determine all others (as was the case in Fig. 4.1 and Fig. 5.3). 
Instead we will see diagrams where three or more states determine all others (as was 
the case in Fig. 5.1b). In these cases calculations will rely more on the analytical 
method of solving the diagrams utilising loops as hrst introduced in section 4.2.5. 
Calculating all loops in a diagram and maximising the argument over these paths is 
a method which works for all diagrams. However, using the method blindly can lead 
to unecessary work, as there are frequently shortcuts to determining the answer, as 
was seen in the case of Fig 5.3. 

The normalised eighth moment is given by 


. tftr('f) 

In this case Wick’s theorem gives (dropping the indices k) 


(5.2.2) 


n n n n 


n n 





n 



tr(I/8) = 2{tTVVVVVVVV) + %{tiVVVVVVVV) + A{iiVVVVVVVV) (5.2.3) 
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n n 

n ^ 

n n 



+ ^tIVVVVVVVV) + ^tIVVVVVVVV) + ^tiVVVVVVVV) 


+ A{tiVVVVVVVV) 


+ ^tIVVVVVVVV) 





1 1 1 

n r 


n 




+ ^tIVVVVVVVV) + lQ{tiVVVVVVVV) 


+ A{tiVVVVVVVV) 



+ ^tIVVVVVVVV) + ^tIVVVVVVVV) + {tiVVVVVVVV) 


(5.2.4) 

(5.2.5) 

(5.2.6) 

(5.2.7) 


(5.2.8) 


where the prefactors indicate the number of equivalent contributions that can be 
obtained by cyclic permutation of the trace. The prefactor 16 in the second term 
of (5.2.6) also incorporates equivalent contributions obtained by reverting the order 
of H’s. Reassuringly the prefactors sum to 105 = (8 — 1)!!, the overall number of 
possible contractions between eight elements. Writing these in terms of (4.2.6) it can 
be shown that 


tr(V?) = 

^^-^ppuu-^uuww^wwvv^vvqq “ 1 “ ‘^^-^ttqq-^qqpp^puvw-^pwvu 
-|- “^ApQiqAqi^yAo^y^^QAcyjr^p -|- ^AipQ^A^yqQAp^yjqA^yjyf 
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(5.2.9) 


The particle diagrams for each of the terms for r are illustrated in Figures 5.5 - 5.14 
and the process for evaluating them is the same as for the fourth and sixth moments; 

(1) diagrams which are otherwise the same as diagrams for lower order moments but 

with the addition of tails simply gain additional factors of for each tail. 

(2) Some diagrams, although initially appearing distinct, actually “collapse” - taking 
the same value as a known diagram - when the condition is imposed on the bonds of 
the diagram that it’s argument be maximal, ie., when we look only at those instances 
of the diagram which attain the maximal argument of the diagram. This process 
will become clearer when looking at some actual calculations. (3) In the absence of 
tails the argument of the particle diagram will be maximised by using the method of 
particle diagrams. This involves determining all loops in the particle diagram and 
identifying which of these, by containing more single-particle states, increases the 
argument of the particle diagram. These methods do not need to be used in isolation. 
For example loop counting can be nsed initially, followed by the observation that a 
diagram will “collapse”, just as it is perfectly legitimate to observe a diagram collapse 
almost fully, except for the addition of a remaining tail. 

A Single Chain 

For the first term AppuuAuuwwAwwwAwqq of r the particle diagram is simply a chain 
of bonds 


p - u - w - V - q 


(5.2.10) 


so that we first choose m states to determine \p) , then add factors of 
each additional tail, giving 



(5.2.11) 
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Diagrams of Lower Order Terms + Tails 

The next term AuqqAqqppApuvw^pwvu illustrated in Fig 5.5 is identical in form to the 
standard diagram but with the addition of two tails, or alternatively Fig 5.1a with 
the addition of a single tail so that 

A A A 

-^ttqq-^qqpp-^puvw -^pwvu 

Likewise the term AeewwAwvupAupqtA^tqw with particle diagram as in Fig. 5.6 is 
identical to the particle diagram of Fig. 5.1b but with the addition of a single tail, 
giving 


'm\ fl — m + 

J I k 


Q Y' ly* 


(5.2.12) 



Figure 5.5: Particle diagram for the term AuqqAqqppApuvwApw^u which is identical to 
the particle diagram of Fig. 5.1a with the addition of a tail, which is the same as the 
standard diagram illustrated in Fig. 4.5 but with two tails. 


AeewwAnjyupAnpqlAyiqyj ^ 


(I — m + k^ 


k 


k 


' I ^ 
,k k k k m — k; 


— k'^ 

V k J 


(5.2.13) 


The same applies to Aqy,ytAuwptApyuqAccpp- Since this is equivalent to the particle 
diagram of Fig. 5.3 but with the addition of a tail we have simply 


A ^ A ^ A A 

■^qwvt^uwpt^pvuq^ccpp 


^m\ (l — m + k'^ 


I 


k 


k 


^ ry^ ry ry ry> ry> ry> 


(5.2.14) 
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Figure 5.6: The particle diagram for the term Ae^wwAwvupAupqtAvtqw which is identical 

illustrated in Eig 5.1b giving 


to the particle diagram of the term AputqAqyj^tAupwv 


this term the same value but with an additional factor of 
given by (5.2.13). 


) (! result 



Figure 5.7: The particle diagram for AtucpAtpcuAcewvAcvwe, although not an extention 
of a previonsly seen diagram with a tail, is simply the prodnct of two standard 
diagrams with a single shared center node given by the state |c). 


The term A^m3p.^tpcu-^cewv-^ 


cvwe 


Finally, althongh we do not have simply a tail extention of a previonsly calcnlated 
diagram for the term AtucpAtpcuAcewvAcvwe illustrated in Fig 5.7, this is a composition 
of two copies of the standard diagram with each copy sharing the state |c). Choosing 
the state |c) will give a factor of after which the reasoning follows identically for 
the left and right instances as it does for the standard diagram so that the binomial 
terms for each are the same as in (4.2.17), giving 


AtucpAtpcuAcewvAcvwe ~ 






(5.2.15) 
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Figure 5.8: Particle diagram for the term AcewuAuwvtAuvqp^pceq of r which takes the 
form of a cuboid structure. The contributing snm is calculated by determining all 
permitted loops in the particle diagram. Then, maximising the argnment over all 
permitted loops reveals the structure exclusively of those loops which maximise the 
argument of the particle diagram (see Figs 5.9a - 5.9f) 








Figure 5.9: Illustration of a selection of the loops within the particle diagram for 
AcewuA^^yfAiiyqpApQ^q illustrated in Fig 5.8. 
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The Cuboid Diagram 


The term AcewuAuwvtAttvqpApceq of (5.2.9) has a cuboid particle diagram, Eig. 5.8. 
One begines by identifying all the permitted loops in the diagram and ends by 
maximising the argument by tuning the number of single-particle states participating 
in each path. The permitted paths for the cuboid diagram are cutpc, ewvqk, cuweh, 
vtpqv, utpqvwu, cutvwec, pcuweq^, cptvqec, cpqeh, utvwu, cuweqvtpc, cutpqvwec, 
utvweqpcu, ptuwvqecp, ce^, mutt, tvt and However, if we check the loops which 
maximise the argument of the diagram we hnd the contributing loops which, unlike the 
permitted non-contributing loops, can contain a non-zero number of single-particle 
states without decreasing the argument of the diagram from its maximal value m — 4k. 
The contributing loops for the cuboid diagram are 


ewvqh 

ce^ 

pq^ 

vtpqv 

utpqvwu 


cutpc 

cutpqvwec 

wuih 

utvwu 

cuweh 

cutvwec 


(see 

Fig. 

5.9a) 

(see 

Fig. 

5.9f) 

(see 

Fig. 

5.9b) 

(see 

Fig. 

5.9c) 

(see 

Fig. 

5.9d) 

(see 

Fig. 

5.9e) 


Moreover, maximising the argument over the contributing loops in analogy to the 
process introduced in section 4.2.5 yields the expression 


AcewuAuwvtAttyqpApQgq ~ 


I 


cm kka^k — ak — f3 



The term A^pcgAg^q^Apu^qA 


uwvt 


We next look at the particle diagram of the term AtpceAevqcApuwqAuwvt of (5.2.9), 
illustrated in Fig. 5.10. To leading order it is the standard diagram in disguise, with 
the addition of 2k degrees of freedom which require explaining. In fact the overlapping 
of the bond t -e with e-c and the overlapping of the bond e-c with c- p 
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on the left sub-diagram must necessarily be minimised in order to maximise the 
argument of the combined diagram. The same applies to the corresponding states 
in the right sub-diagram of Fig. 5.10. The minimum overlap is of course m — 2/c as 
before, giving 

e - V C e - 1 (5.2.17) 

o- -q C c- p (5.2.18) 

e-c C e-c (5.2.19) 

The additional coincidence in m — 2k states between say |f) and \p) which already 
share k single-particle states through the bond t --p brings the total overlap between 

/ 7 \ 2 

them to m — k. This gives an additional factor one factor each for the left 

and right sub-diagrams which together form Fig. 5.10. One also gains the double 

factor from the choice of single-particles contained in e-c fl e-c and 

correspondingly u - w fl u - w. In this way we find 


^tpce-^evqc-^puwq-^uwvt ^ 


^ I ^ 
,s k k k k, 


m 


kV (I 


2k 

k 


(5.2.20) 


This is an instance of a particle diagram “collapsing” partially to a diagram seen 
before in the calculation of a lower moment, in this case the standard diagram of Fig 
4.5 which has been a presence in calculations from the fourth moment onwards. 


The term ApetqAq^^vAgu^cA 


cwvp 


The calculation for the term ApetqAqtuv^euwc^cwvp of (5.2.9) with diagram shown in 
Fig. 5.11 will be calculated using the method of particle diagrams by utilising loop 
summation. Explicitly, the number of single-particle states contained in each loop is 
optimised in order to maximise the argument. The contributing loops are those which 
do not necessarily have to contain zero elements in order to maximise the argument 
given by the sum arg = where rii is the number of elements in loop i. The loops 
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®-© . © . ® - © 


®- © . ® (£> . @ - ® 

Eigure 5.10: The AtpceAevqcApuwqAuwvt which to leading order is the same as the 
standard diagram, bnt with an additional k degrees of freedom. This is analogons 
but not exactly the same as the case where we have a known diagram with an extra 
tail. Eor k < m — k the overlap of the bonds linking the states |e), |c) and |m), \w) 
respectively implies that the states |t) and \p) overlap by an additional s = m — 2k 
states, bringing their total overlap to m — k and likewise for the states |t) and |n). 

@ - ® © - © - @ 

I II I I II I 

I II I I II I 

I II I I II I 

I II I I II I 

© © - © © @ © 

Eigure 5.11: Illustrated is the particle diagram of ApetqAqtuvAeuwcAcwvp which has 
some similarities to the standard diagram. Eor the purpose of evaluation however, 
it is not similar enough. To calculate the leading order summation for the terms 
satisfying this diagram one needs to use the method of particle diagrams, availing of 
loop summation. 


on the particle diagram of the term ApetqAqtuvAeuwcAcwvp are 


a = pv~ 

(3 = pe5 

7 = 


^ = qvpeut 


9 = qpvuet 
A = cpeuvib 
5 = uvih V = cepvuiJb 


e 

V 

00 

p 


= qtuv 
= qtep 
= cwv^ 
= cwuh 


71 

a 

T 

P 


pqvwutec 

pvuk 
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Because the number of single-particle states contained in all paths which pass through 

a bond-must sum to k and similarly the number of single-particle states contained 

in all paths which pass through a bond-must sum to m — k one can immediately 

read off the following conservation equations 

Ua + + ng + n,, + Up = k 

71/3 + + nx + n,, + Up = k 

+ ng + nx + Up = k 

ns + ng + nx + n„ + np = k ( 5 . 2 . 21 ) 

= k 

Up + ng + n„ = k 

Uaj + nx + Ur = k 

Up + Uiy + Ur = k ( 5 . 2 . 22 ) 

Uf: + + n^j + Utt + Ua = m — k 

He + + Up + tIt, + = m — k 

+ nx + Ut, + = m - k 

+ Up + riy + Ut, + ns = m — k ( 5 . 2 . 23 ) 

Up + ng + Uuj + Ut, + Ua = m - k 

Up + + nx + n-j, + np = m - k 

Up + Up + Uu + Ut, + np = m - k 

Up + ng + Up + + Up = m — k ( 5 . 2 . 24 ) 

Eor a particular loop to contribute the argument must additionally satisfy 

arg = ^Ui = Ua + np + Up + ns + + Up + + Up + + ng + nx + n,y + + Up + + Ur 



5.2. THE EIGHTH MOMENT 


113 


= m + Ak (5.2.25) 

since m+Ak is the maximal argument that this diagram can attain and therefore the rii 
satisfying this equality are the only terms of interested in the limit / —>■ oo. Imposing 
this contraint on the conservation equations yields Ua = rij^ = = k and 

= ne = nx = n,y = Hp = 0. Moreover 


Tl^ ^rf) 


T^tj- Tl'j^ 

Ua = Ua 
Ur = fir- 


Inserting these back into the conservation equations gives 


= k — Her 
ni_j = k — Hr 

+ Uuj + n-jr = m — 2k 


(5.2.26) 

(5.2.27) 

(5.2.28) 


so that in the limit I —)■ cx) the surviving term reads 


-^petq-^qtuv-^euwc-^cwvp ~ ^ ) 


= E 


^ I ^ 

^Ua ng ns n„ Ur n^ 

I 


y,n^ 


,k k k k n„ Ur k — n„ k — n„ k — Ur k — Ut- m — 4k + n„ + Uj 


(5.2.29) 


Dividing (5.2.29) by the normalisation term (^tr(I4^)^ and summing out r one 
extraordinarily attains the following Hahn polynomial 


T 4 A A A 
j. jY ^petg^gtuv^euwc^cwvp 

N—^od 


(if(i^)) (:) 


(y 9 ^ /m — k — I'm — 2k^ ^ k 


k 


Un 


Un 


(5.2.30) 


This is unusual because no previous term has explicitly required a Hahn polynomial 
in order to be expressed properly. Additionally, this is an interesting development 
because it hints that it may be possible to express all previous terms as Hahn 
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polynomials and perhaps even find a set of tools with which to calculate the moments 
using a “Hahn Method” as apposed to using the method of particle diagrams which 
this thesis sets out to present and explain. 


5.2.1 Hahn Polynomials 

The author’s proof of the following lemma will briefly be highlighted to express (5.2.29) 
as a quotient of binomials in m and k. The proof uses the formalism of the class of 
Hahn polynomials for which definitions can be found in the compendium [KRLSIO]. 

Lemma. 


(m - k - a\ f k\ 
k )[a) 


m 


2k'^ 


a 



(5.2.31) 


Proof: 

Writing the l.h.s as an hypergeometric function 


E 


— k — a\ (k\ (m — 2k'' 


k j \aj \ a ) \ k ^ 

and recalling the definition of a Hahn polynomial 


m — k\ I —k,2k — m,2k — m 

3-^2 I _ ; 1 


l,k — m 


n t n Ar\ . -n,n + a +/S + 1 ,-x 

Qn{x-,a,^,N) := 3F2 I ^ ;1 

a + 1, —N 


(5.2.32) 


(5.2.33) 


we have n = k,a = 0, fi = k — m — 1 and x = m — 2k. To express this as a series it 
should then be noted that 


iFi 


—X 

o T 1 



iFi 


x-N \ 

/9 + 1 ’7 


E 


i-N)n 
(/3 + l)nn! 


Qn{x;a,/3,N)F. 


(5.2.34) 
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(see [KRLSIO] chapter 9). Furthermore, by dehnition 


a 

b 


where 


[a := 


r(a + n) 

r(a) 

and the symmetric extention of this is 


\ - (a)„r 

) ■ h'o ib)n n\ 

(5.2.35) 

(ft “h 1)(^ “h 2) . . . (ft “1- Tl) 

(5.2.36) 

l)"'{x -n + 1)„. 

(5.2.37) 


Using the dehnition it then follows straightforwardly that 


.ni (8.2.38) 


nUo r(m — 2fc + 1 — n)n\ n\ 


n=0 


n n\ 


and similarly 


iF, 


-k \ ^ T(k + l)T(m-k-n+l) f‘ ^ Q t" 
-{m-k)j ^QT{k-n + l)T{m-k + l)n\ n\ 


so that the complete series can be written as 
2k — m 


00 


1 ’ (V) 


Finally, comparing the coefficients of gives 


(5.2.40) 


1 ^ /fc\2/m-2fc\ 

[ P ) 


(5.2.41) 


which completes the proof. 


□ 
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5.2.2 Remaining Diagrams 

There are just three remaining terms left in the calculation of (5.2.9) for which we have 
already fully covered the techniques required; (1) comparison to diagrams representing 
lower order terms, (2) the “collapse” of diagrams to others of familiar form such as in 
the case of Fig 5.10 which collapses to take the form of the standard diagram with a 
tail when we just look at those components of the diagram which attain the maximal 
argument, and hnally (3) the method of particle diagrams in which we maximise the 
argument ( = degrees of freedom) of the diagram by identifying all permitted loops 
in the diagram and from these finding all contributing loops. The path snmmation 
method can be used in every case, even those where (1) and (2) do not apply. 


The term Au^vtA^qvtApueqAp^ec 


Identifying paths and maximising the argnment for Auwvt^cqvt^pueqApwec of (5.2.9) 
with corresponding diagram Fig. 5.12 in the same manner as described previously we 
yields the following paths 


a = utvwe 
(3 = tcpqv 
7 = ueih 
6 = 

e = utce 
r] = wvqk 
u = qvwj) 
fi = ctup 
^ = eutcpqvw 
p = eutvwpqc 


a = ewvtupcq 
T = pqvtceuw 
(j) = pctvqewu 
9 = ectvq 
X = putvw 
V = ec^ 

TT = pw^ 

'ip = euwpqh 
K = pcqewii 
X = ^ 


and the following loop conservation eqnations 


Ua + nx + Ucr + Up + = k 
ng + ne + Ut- + n^ + = k 


(5.2.42) 
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+ ne + + Up + Ur + = k 

Up + ne + = k 

Uuj + nx + Ut, + Up + Ur + = k 

Up + nx + + Ua + = k 

Uy + Ut, + Ur + = k 

ns + n„ + Up + = k ( 5 . 2 . 43 ) 

Ua + n^ + n^ + + Up + n-r + = m — k 

Ua + + Up + + Ua + + Uf^ = m - k 

Ua + Ue + Up + + nx + Up + Ua = m - k 

Ua + Up + Uaj + + nx + Up + Ua = m - k ( 5 . 2 . 44 ) 

ni 3 + + Up + + ne + Ur + = m - k 

nis + ns + Up + = m - k 

np + ns + n^ + n^ + np + nr + n^ = m - k 

np + np + n^i + n^ + ne + nr + n^ = m — k. ( 5 . 2 . 45 ) 

Solving for 

arg = ^ni = na + njs + np + ns + ne + np + n^ + np + ni^ + ne + nx + ny 

i 

+ n^ + np + n^ + nr + n^ + n^ + n,^ + n^ = m + Ak ( 5 . 2 . 46 ) 

gives the following restrictions on the number of single-particles in each of these loops 

ne = nx = ny = nr = np = n„ = nr = n^ = n^ = n^ = 0 ( 5 . 2 . 47 ) 

np = ns = n^ = np = n^ = np = k ( 5 . 2 . 48 ) 
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Eigure 5.13: The particle diagram for Auwce^epqc^tqvu^twvp- By acknowledging the 
overlaps between the bonds of the left sub-diagram we can conclude that the complete 
diagram is identical to Fig 5.3 which we have already calculated for the sixth moment, 
with an additional factor degrees of freedom between |c) 

and |e). 


with Ua = Up, as well as the following identities 


= k — ria 
= m — 3k — ric 


(5.2.49) 

(5.2.50) 


Hence 


AuwvtAcqvtApueqApyj^Q ~ ^ ) 




) 

k k k k k Ha np n^J 
' I 

,k k k k k k Un rin k — Tin m — 3k — rir, 


(5.2.51) 
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(d) The loop tqpvuecwt 


(e) The loop tqcwt 


(f) The loop wvqcil) 



(j) The loop twut 


(k) The loop qvp}^ 


( 1 ) 


The loop twceut 
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© 

(a) The loop cqvpeh 




(d) The loop tqcept 







@ © ® 

® [] ® 
® ® ® 


(g) The loop wvqptuw (h) The loop ce^ 

Eigure 5.13: Illustration of all the loops in the particle diagram for the term 
Auwce®pqcAqvuAwvp represented in Fig 5.13. 
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giving 


lim 

N^oo 


N 


■uwvt -^cqvt -^pueq -^pwec 



(5.2.52) 


Th0 torin Auwce^epqcAtqvuAtwvp 


The penultimate term AuwceAepqcAtqvuAwvp 
Fig. 5.13 which has the particle loops 


of (5.2.9) yields the particle diagram of 


a = wtqcepvu 
(3 = wceutpqv 
7 = wutpecqv 
6 = wceuvpqt 
e = wtqc 
p = wvqh 
u = tuep 
fi = ueph 
^ = twcqvpeu 
9 = wtu 


A = qv^ 
p = wvtt 
a = 

V = wtuec 
71 = cqvph 
T = tqcep 
(p = wceul 
pj = wtqpvu 
K = wvqptu 



and the following particle conservation equations 


Uq 

, + ns 

+ n^ 

+ 

rig 

+ 

ne 

+ 

n^ 

+ 

np - 

= k 

(5.2.53) 

Ua 

+ ns 

+ np 

+ 

ng 

+ 

nx 

+ 

nn 

+ 

np - 

= k 

(5.2.54) 

Ua 

+ n^ 

+ n^ 

+ 

np 

+ 

ng 

+ 

n-n 

■ + 

n-r = 

= k 

(5.2.55) 

Ua 

+ n^ 

+ n^ 

+ 

np 

+ 

ng 

+ 

n-n 


n-r = 

= k 

(5.2.56) 


Ua 

+ ns 

+ 

n. 

+ 

n-r 

+ 

Tl^ 

+ 

n„ = 

= k 

(5.2.57) 


na 

+ ns 

+ 

np 

+ 

np 

+ 


+ 

np - 

= k 

(5.2.58) 


Ua 

+ n^ 

+ 

no 

+ 

np 

+ 


, + 

n^ = 

= k 

(5.2.59) 



na 

+ 

n^ 

+ 


+ 

Ur 

+ 

nx = 

= k 

(5.2.60) 
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rip + + ne + = m — k ( 5 . 2 . 61 ) 

rip + rip + n\ + Ht, + rif^ = m — k ( 5 . 2 . 62 ) 

np + ns + + rir^ + np + riy + rip = m — k ( 5 . 2 . 63 ) 

np + ns + Uui + Ufj, + np + Uy + np = m - k ( 5 . 2 . 64 ) 

np + n^ + n^^ + np + np + n^ = m — k ( 5 . 2 . 65 ) 

np + n^ + n^ + n„ + rir + nf^ = m — k ( 5 . 2 . 66 ) 

np + ns + nx + n„ + np + n^^ = m - k ( 5 . 2 . 67 ) 

np + ns + ny + np + n-^ = m — k. ( 5 . 2 . 68 ) 


Solving for 


arg = ^ni 

i 

= na + np + n^ + ns + n^ + np + n^ + np + np + ne + nx + n,y 


+ n.„ + np + n„ + n^ + np + np + n^ + n^ 

= m + 4:k ( 5 . 2 . 69 ) 

gives the following identities 

np + n^ = rn — 3k ( 5 . 2 . 70 ) 

ne + niy = k ( 5 . 2 . 71 ) 

n0 + n^ = np ( 5 . 2 . 72 ) 

np + nf, = n„ ( 5 . 2 . 73 ) 

as well as the restrictions na = n-y = ns = n^ = n^ = np = n-j^ = nr = np = 0 


and np = n^ = nx = n„ = n^y^ = k on the number of single-particles in each of 
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these loops. In addition one finds 


710 

= m — 3k — 


He 

= n^-n^ = k-rip-n^ 


riy 

= rip+ 71,, 



= k — rip 

(5.2.74) 

Up 

= rip 


rin 

= 

(5.2.75) 


so that the final expression takes the value 

I 

kkkkkKm — 3k — upp + nk — p 


A A A, 

^uwce-^epqc^tqvu^twvp 




which then gives 


lim 

N^OO 


2-4 


uwce^epqc-^tqvu^twvp 





(5.2.77) 


The final term Auvqt^^twvc^^cewpApueq 

Eor the evaluation of the hnal term AuvqtAtwvcAcewpApueq of (5.2.9) one can use the 
same method as outlined in the evaluation of ApyuqAg^^tAtpwu (Fig- 5.3). That is, 
note that the states \u) and |u) together determine all other states in the diagram, 
Fig 5.14. This means that to maximise the argument of the diagram we minimise the 
overlap between \u) and |u). As before this is done by identifying as many loops as 
possible within the diagram from \u) to |m) without passing through |u) and from 
|u) to |u) without passing through |m). The single-particle states which cannot be 
incorporated into these loops must necessarily be in the set that is shared by both \u) 
and |u). And whereas for the case of Eig 5.3 each loop is composed of three bonds, 
for this analogous term of the eighth moment each loop consists of | = 4 bonds. The 
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Eigure 5.14: The particle diagram for Auvqt^twvc^cewp^pueq which takes the form 
(skkklkkkk) with s = m-4k. 


loops passing through the state \u) without passing through |u) are 


Pi = utcpu 
P 2 = utweu 
P 3 = uewvtt 
P 4 = utcpweu 


(5.2.78) 

(5.2.79) 

(5.2.80) 

(5.2.81) 


and similarly the loops passing through the state |m) without passing through |u) are 


P5 = qpcv}j[ 

(5.2.82) 

Pe = qewv^ 

(5.2.83) 

P7 = qpctq 

(5.2.84) 

Ps = qvwecp^. 

(5.2.85) 


The single-particle states which are not an element of each of the bonds in a loop 
are necessarily an element of all the bonds in a loop through |m) and |n) and there 
are s := m — 4k of these. The final expression then reads 

AuvqtAtwvcAcewpApueq 


(5.2.86) 
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Bringing all terms (5.2.11 - 5.2.16), (5.2.20), (5.2.30), (5.2.51), (5.2.76) and (5.2.86) 
together gives 



Illustrated in Eig 5.15 the eighth moment takes the same form as the fourth and 
sixth moments, transitioning from a gaussian value (2n — 1)!! = 105 at k = 0 to a 
semi-circular moment for all fc > y. Unlike a the lower moments the eighth moment 
contains a Hahn term. Fortunately the method of particle diagrams can be used to 
yield all such terms, and as we have seen the procedure is the same, independent 
of how complex a particular diagram is. These results and the method of particle 
diagrams will be discussed in greater detail in the next chapter. 
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Figure 5.15: The eighth moment r of the level density against k/m for m = 6,8, 10,12 
exhibiting the same properties as the fourth and sixth; a semi-circular domain, r = 14 
for k/m > ^ converging to a gaussian value, r = 105 for k = 0. Higher values of m 
give faster convergence to the semi-circular moment. 





Chapter 6 


Conclusions 


6.1 Moments 

This thesis has focused on the problem of characterising the statistics of unified RMT, 
more commonly known as embedded many-body random matrix theory. As noted 
in chapter 1 this is a relatively yonng area of research, with few practitioners and a 
modest body of literatnre. The approach established here has involved the method 
of particle diagrams, a mathematical toolkit invented by the author, with which we 
have investigated the calculation of moments of the level density p{E) of the random 
potential derived in chapter 2 given by 

Hk = '^Vjia]ai. ( 6 . 1 . 1 ) 

ji 

These potentials are characterised by the parameter set {k, m, I, S'} where k is the 
order of the potential, m is the number of particles in the experiment, I is the number 
of single-particle energy levels and S represents any additional symmetry conditions 
placed on the potential H^. Together these determine the state space of the potential, 
and thereby the state space of unified random matrix theory. This is a superspace of 
canonical random matrix theory. Canonical RMT places a single random variable into 
each matrix cell and coincides with the point {m, m, oo, S}, with S usually enforcing 
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orthogonal (/S = 1), unitary /? = 2 or symplectio (/I = 4) symmetry on Hk- Here we 
have made calculations of the fourth, sixth and eighth moments for all fc < m for 
the embedded CUE (see section 4.2). It was noted in (4.1.5) that the moments of a 
Gaussian are given by 





= (2n - 


1 )!! 


( 6 . 1 . 2 ) 


and the moments of a semi-circle, those characterising Wigner’s semi-circle law, are 
given by 


— x‘^ dx = Cn 

Tir^ J-r 

where the Catalan numbers (3.2.57) are 


(6.1.3) 


a 


1 

n -I- 1 



(6.1.4) 


Either through direct calculation or by using the method of particle diagrams we 
have seen in section 4.2 that the fourth moment of the level density is 


lim K 

N—^OO 


A 4 A 
2 lim jv 

N —^■CXD 1 A 

V A 





(6.1.5) 


which is illustrated in Fig 4.2. The hgure reveals a transition from a Gaussian fourth 
moment {2n — 1)!! = 3 at fc = 0 to a semi-circular moment Cn = 2 for all /c > y. 
This includes k = m which we know already from Wigner’s semi-circle law should 
yield k = 2. Hence there is a portion of the domain in the unihed theory (y,m] in 
which k varies (and therefore the action of the potential Hk changes) but the fourth 
moment remains fixed to a semi-circular value. 

It is shown in chapter 5 that the sixth moment is given by the combinatorial sum 

fm—k\ fm—2k\ fm—k\ /m— 

lim h = 5+ ^ ^ +6VM- + 3 Vx 2 • (6-1-6) 

N—^oo /m\^ im\^ 
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This was illustrated previously in Fig 5.4. Again we can observe a gradual transition 
from a Gaussian valued moment (2n — 1)!! = 15 at /c = 0 to a semi-circular moment 
Cn = ^ not just for the canonical case k = m, but for all A; > y. 

The eighth moment, also derived in chapter 5, was shown to be equal to the 
rambling expression 



which has been illustrated in Fig 5.15. From the figure we can immediately note the 
way in which the value of the moment transitions from a gaussian value (2n—1)!! = 105 
on the left-most edge of the phase space /c = 0 to a semi-circular moment = 14 for 
all values y < fc < m where despite the potential Hk being fundamentally different, 
the eighth moment remains constant. 

6.2 Internal Structure of the Moments 

The results shown above suggest that the unihed RMT state space for the eGUE 
consists of two distinct domains in the order k. Firstly, each of the moments reveal 
a canonical domain dehned by y < fc < m (see Fig 6.1). In this region of the 
phase space the action of the potential is changing, but the moments of the level 
density remain fixed to semi-circular values. Hence, we may expect that in this region 
Wigner’s semi-circle law applies for all values of k. 

Secondly, again with reference to Fig 6.1, there is a critical domain dehned by 
0 < A; < y. Here the value of k is crucial for determining the value of the moments, 
which vary for all k within this region of the state space. 
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Figure 6.1: Illustration of the general form of all the moments, showing a transition 
from a Gaussian value to a semi-circular value for all k in the canonical domain, and 
some additional structure in the critical domain at the discrete points k = 


Finally, the critical domain for each of the moments contains an internal strnc- 
tnre. For the 2n’th moment this additional detail is dehned by the discrete points 
y, y,..., ™ at which the combinatorial expressions (6. 1.5-6. 1.7) contain terms which 
become zero for all valnes of k thereafter. For example the qnotient 
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giving the contribution of the standard diagram to the fourth moment (6.1.5) will be 
zero for a\\ k > m — k while the quotient 


( m—k\ / m—2k\ 

k ){ k ) 



( 6 . 2 . 2 ) 


from (6.1.6) will be zero for all k > m — 2k. The gradual removal of these terms at the 
discrete values of fc = y,..., ™ is what produces the decline in the moments from 
gaussian values to semi-circular values and explaining the physical ramifications of 
this (if any) is an open question which would make for some very interesting reading. 
Related to this, but offering perhaps less detail, is a physical explanation for why the 
random potential, being different for each value of k, should nonetheless have the 
same level density for all k in the canonical domain k > ^. 

A fundamental component of the internal structure of the critical domain is the 
existence, in each of the moments, of a single diagram which takes the value 


^ / 

,s k k ... k 


2n times 


(6.2.3) 


where for the 2n’th moment s := m — nk and the sequence of fc’s repeats 2n times so 
that the sum is s -f J22n k = m + nk. This is the argument attained by any diagram 
contributing to the limit value of the 2n’th moment as / —)■ cx). For the fourth moment 
the diagram which takes the form of (6.2.3) is the standard diagram. Fig 4.5, while 
for the sixth and eighth moments it is Figures 5.3 and 5.14 respectively. The fact that 
these diagrams take the form of (6.2.3) is due to the existence in each case of just 
two m-body states which together determine all other single-particle labels present in 
the diagram. 
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6.3 The Dilute Limit k m <^l 


Another result of the analysis is the conclusion that in the limit case k/m ^ 0 
the moments are those of a Gaussian distribution. This can be seen in equations 
(6. 1.5-6. 1.7) by taking the limit m —)■ oo and applying Stirling’s approximation to 
the resulting binomials, yielding the moments for the special case k <^m <^l which 
is also known as the dilute limit. In this limit each binomial expression contributes 
just a factor 1, so that the moments become equal to the sum of the coefficients. 
Since the coefficients sum to the number of unique pairwise partitions of a set of 2n 
objects the result is always (2n — 1)!!, the Gaussian moments. The fact that each 
combinatorial term gives a contribution of unity can also be seen by noting that the 
number of combinatorial expressions in m in the nu m erator is always equal to the 
number of combinatorial terms in m in the denominator. For example the component 

( m—k\ {m—2k\ / m—3k 
k )\ k )\ k 

& 



of the eighth moment (6.1.7) has three binomial terms in nr 

forming the numerator and three binomial terms (jfj in m forming 

the denominator. Likewise, the Hahn expression 


2 


fc °) fm-k-a\ 

{-y ^ ) 




(6.3.2) 


has the three factors ^ in m forming the numerator and the three 

factors in the denominator. 


6.4 Bosonic states 

Using the method of particle diagrams it is possible to make a hnal remark regarding 
the statistics of the hermitian eGUE potential (4.2.4) given a bosonic state space. 
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From section 2.1.3 it is known that bosonic many-particle states may contain multiple 
copies of the same single-particle states. Retaining the same potential (4.2.4), but now 
assuming a bosonic state space, it is possible to use what has been learned already 
for the eGUE embedded in a fermionic state space to make far-reaching conclusions 
about the bosonic case. Namely, it is straightforward to see that in the limit I —)■ oo 
the results for all moments agree with fermionic systems. Intuitively, this happens 
because in this limit contributions arise only for those choices of many-particle states 
which maximise the number of participating single-particle states. This means that 
in the bosonic case multiple occupancy of the same single-particle states is penalised, 
such that for the terms which survive in the limit I —)■ oo there is no difference from 
the fermionic case. 

Eor a formal derivation consider bosonic m-particle states containing repeats of 
z < m unique single-particle states. There are ways to select the participating 
single-particle states, and ways to select their multiplicities to obtain altogether 

particles. As stated before in section 4.1.3 the number of many-particle states is 
thus 


^ = E 






(6.4.1) 


where the overall sum as well as the summand z = m have the argument m whereas 
all other summands have lower arguments. As a consequence the asymptotic form of 
N coincides with the fermionic case where N = 

The same logic applies when evaluating the particle diagrams, where the dashed 

bonds-contain only unique sets of k single-particle state labels whereas in 

contrast to the fermionic case the solid bonds-contain m — k single-particle state 

labels with repeats now permitted. In the fermionic case, the contributions to the 
particle diagrams with maximal argument m nk can always be factorized into a 
term counting the number of ways in which the m + nk participating states 

can be chosen, and further /-independent factors counting the number of ways in 
which these m + nk states can be distributed among different sets while obeying the 
conditions implied by the diagram. In the bosonic case we instead have to select 
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1 < z <m + nk states with multiplicities summing to m + nk. The contribution of 
each diagram turns into a sum over and for each ^ a single /-dependent term 
is obtained, together with further /-independent hnite factors counting the number 
of ways to distribute these states among sets containing m particles. The latter 
choices have to be compatible both with the restrictions implied by the diagram and 
the requirement on the multiplicities. However, the only summand attaining the 
maximal argument will be the last, for which z = m + nk. For this summand all 
multiplicities are 1 meaning that we only have to consider bosonic states composed 
entirely of distinct single-particle states. Since these are mathematically equivalent to 
fermionic states it can be concluded without any further calculations that all bosonic 
and fermionic moments are equal in the limit / —>■ cx). 

6.5 The Method of Particle Diagrams 

The focus of this thesis has been the study of quantum many-body potentials taking 
the second quantised form of (6.1.1). This involved the study of a unihed phase 
space of random matrix theories composed of the parameters k,m,l which are the 
order of interaction, number of particles and number of available states respectively. 
These, in addition to the symmetry constraints on Hk, determine each point in 
the phase space. The canonical form of random matrix theory corresponds to the 
single point {m,m,oo}. In their seminal paper Benet, Rupp and Weidenmiiller 
were able to calculate the fourth moment of the level density of the eGUE in the 
domain {0 < k < m,m, cxo} [BRWOla]. In so doing they also illustrated the difficulty 
involved in studying the more general embedded RMT state space, suggesting that 
new methods would be needed beyond those which are specialised for the canonical 
k = m case. The author’s response to this was the creation of the method of 
particle diagrams which has been illustrated in the previous chapters to calculate the 
fourth, sixth and eighth moments of the embedded GUE. This is a modest success 
in itself, since the sixth and eighth moments were not known before in the critical 
domain k < ff. Moreover the method is unmatched in it’s simplicity, involving no 
superalgebra and no complex mathematics beyond basic combinatorics. The basic 
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ingredients are as follows. 


1 . 


Arguments (order of magnitude). The first ingredient of the method of particle 
diagrams is the quantity called the argument which is the power in / of a binomial 
expression and can be shown using Stirling’s formula limn^oo'^^- = A/27rn(^” j 
to be given by 


arg 




(6.5.1) 


Arguments are a simplifying feature of the method of particle diagrams, working 
to seperate those terms which will not contribute to the limit value of the 
moments as / —)■ oo but simultaneously leaving enough mathematical structure 
behind so that the limiting result can be calculated. Arguments also provide 
the constraints neccessary to optimise loops. 


2. Particle Diagrams (A graph structure to represent factors of A^^p^). The 
particle diagrams, although themselves not essential to the calculations, aid us 
in visualising “what is going on” as we proceed to calculate the arguments and 
identify loops in the diagram for each term of the higher moments. Particle 
diagrams form the bridge between arguments and loops. They also provide a 
way of identifying visually, as opposed to analytically, all the loops in a given 
term. Particle diagrams found in calculations for lower moments re-appear in 
calculations for higher moments so that calculating successive moments makes 
later work easier. It also means that the standard diagram. Fig 6.2, is found in 
calculations for all the higher moments, even though it’s maiden appearance is 
in the calculation of the fourth moment. 


3. Loops (Intersections between states). Loops represent intersections between 
states. These are the single-particle state labels shared between the many-body 
states (which are sets of m state labels). Loops also represent the sets which 
must be summed over in order to calculate the trace which yields the moments 




(6.5.2) 
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Figure 6.2: The particle diagram of the term is the standard diagram. 

Since diagrams from lower order moments appear again in calculations of higher 
moments the standard diagram appears in calculations of all the moments. The 
standard diagram is also an example of a diagram whose value takes the form of (6.2.3) 
since the states |/r) and |p) together determine the single-particle states contained in 
|cr) and |i/). 


Finally and most importantly, loops determine the argument of a diagram, so 
that the sum of loops in a diagram must be maximised. 

Together these three ingredients form the method of particle diagrams. Although the 
method has been used here to calculate moments for the eCUE it can be used to 
calculate statistics for other symmetry classes as well. The most interesting of these 
is the embedded gaussian symplectic ensemble (eCSE) about which virtually nothing 
is currently known. The author has also calculated the fourth and sixth moments for 
the case where the potential satishes the eCOE symmetry of (4.1.21) and moreover is 
embedded in a fermionic state space. In this instance it is seen that, conditional on 
taking the limit / —)■ cx), these moments are identical for the eGUE and eGOE. 

To end, it should be noted that although particle diagrams are used here to 
calculate moments for all 0 < fc < m, the results coincide with Wigner’s semi-circle 
law as expected for the canonical case A: = m as well. In fact, for k = m the dashed 

bonds-must necessarily contain k = m elements, whereas the solid bonds- 

contain m — k = 0 elements. This in turn results in diagrams which are equivalent 
to those shown in section 3.2.3 where cycles are represented as partitions of regular 
polygons, or equivalently Dyck paths. 
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6.6 Related Literature and Wider Context 

6.6.1 Many-Body Potentials Revisited 

The unification of random matrix theory offers the possibility that the landscape 
of statistical qnantities of random matrices can be stndied en masse, with a single 
eqnation for each qnantity. Each statistical qnantity found in canonical random 
matrix theory would then be just a special case {k = m) of an equation determined 
by the parameters k, m and I as well as the symmetry conditions placed on the 
matrix elements. There is another paradigm, however, which is larger still than 
the embedded state space and which amalgamates and connects a vast array of 
physical and mathematical models. This is the paradigm of many-body potentials. 
The set of many-body potentials includes random matrix theory, both its canonical 
and unified form (6.1.1), as well as conntless other models. These models do not 
share the same hamiltonian but are unified by the similarities in their mathematical 
strncture. Together they explain a diverse array of disparate phenomena from the 
study of quantum spin hypergraphs to the application of neural networks implemented 
algorithmically by compnters (essentially “artificial brains”) which underly most 
modern online services from Google, Amazon, Netflix and countless other technology- 
driven enterprises. The fnndamental form nnderlying all many-body potentials is 
given by the equation 

( 6 . 6 . 1 ) 

kik2 

where |fci| = |fc 2 | = k are fc-tnples. Here and 3^2 are any fc-tuple operators. 
They can be complex and mnltidimensional so, for example, they can be matrices. 
The coefficients v{ki,k 2 ) are numbers, possibly complex and if they are random 
variables one refers to (6.6.1) as a random many-body problem. The dehning feature 
of many-body potentials is the fact that the operators and 5^2 are fc-tnple 
operators. 
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m-tuples 





^-tuples *0 • • 

• • • 


Figure 6.3: Illustration of a many-body problem, given by (6.6.1). The hamiltonian is a 
sum over fc-tuples in an m-tuple with k < m. The m-tuples usually represent quantum 
many-body states, and the fc-tuple operators with coefficient v{ki, ^ 2 ) give 

the interaction terms between fc-tuples of single-particle states within the compound 
m-body state. 


6.6.2 Other Unified RMT Results 

The embedded form of RMT (6.1.1) introduced by Mon and French [MF75] coincides 
with the general form (6.6.1) when the fc-tuple operators 

Aki = aj 

Bk, = ai ( 6 . 6 . 2 ) 

and v{ki,k2) = Vji where Vji is a Ganssian random variable and as before j = 
(ji, • • • ,jk)- These ensembles have been studied before by [MF75,FW70,BRW01a, 
Kot03, N014, Sre02] among others. For the limit case k/m —)■ cx) it was fonnd 
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by [MF75] that the distribution is a Gaussian. This is corroborated by the results as 
discussed in section 6.3 above. Additionally, before the method of particle diagrams 
was used to calculate the same results Benet, Rupp and Weidenmuller (BRW) had 
already shown using their eigenvector expansion method that the fourth moment of 
the eGUE is given by (6.1.5). They also showed using a mixture of supersymmetry and 
eigenvector expansions that the moments are semi-circular in the canonical domain 
fc > y [BRWOla]. The method of particle diagrams, which is the latest addition to 
the family of embedded many-body methods, reprodnces the fonrth moment of the 
eGUE as shown by BRW and adds the sixth and eighth moments for all 0 < A; < m. It 
also reveals the internal structure of the moments in the critical domain as discussed 
in section 6.2. 


6.6.3 Quantum Spin Graphs 

An area of current research which yields fascinating similarities to the embedded 
many-body state space is the study of quantum spin graphs. These are many-body 
potentials taking the form 


T, E 

(jl,- 5fc)er„ ai,--- ,afc=l 


a. 


0 - 1 , 


,Ok{jl 




(“l) 

ii 


' a 


(Ok) 


3k 


(6.6.3) 


Gomparing this with (6.6.1) shows that the fc-tnple operator is 

/ID _ ^Ul) nU*) 

Aki^k2 <^ji '"(^jk 


(6.6.4) 


and the coefficients are 

u(fci,fc2) = aai,-,ak{3l,-,3k)- (6-6.5) 

In this context k is the number of vertices in every hyperedge of the graph T. Viewing 
(6.6.3) in the context of the general form (6.6.1) shows us that this is a generic 
many-body hamiltonian with Pauli spin tuples acting as the /c-tuple operators. It is 
also known as a quantum spin glass model. If A; = 2 the potential coincides with a 
quantum spin chain model. These hamiltonians, particularly the form with k = 2, are 
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the focus of a growing body of research which includes [ES14,HMH05,GFS12,KLW14]. 
The work of [ES14] is particularly interesting in the context of this thesis, since it was 
proved there that the level density of a random hypergraph (6.6.3) takes the form 
of a Gaussian distribution for k <C \/m, a semi-circle for k 3> ^/m and a compactly 
supported density function for values of k and m satisfying Urrim^oo^^ = A > 0. 
This shows that the statistical features of quantum hypergraphs are strikingly similar 
to the results shown here for the unihed RMT state space. Both show a transition 
from a Gaussian regime to a semi-circular one. While the point of transition in the 
case of the eGUE is fc = y, the results of [ES14] suggest it is at fc ~ ^/m for the 
quantum spin hypergraph potential. Finally, whereas for the eGUE one takes I —)■ cxo, 
for spin graphs there are only two available states, up and down, so that I = 2. Hence, 
although the embedded RMT potential (6.1.1) and the quantum spin graph potential 
(6.6.3) are clearly different, they also share many statistical features and both models 
take the form of a generic many-body potential (6.6.1). Exploring their statistical 
similarities further could yield some remarkable insights into many-body potentials. 

6.6.4 Neural Networks 

Many-body potentials of the form (6.6.1) are also present in research into artihcial 
intelligence, notably machine learning neural network models which attempt to 
simulate some basic functionality of a generic brain. These machine learning networks 
are used by the vast majority of modern technology companies such as Google (search), 
Amazon (online shopping), Netflix (movie streaming) and IBM (most famously in 
it’s Watson machine which is now the world champion of the popular quiz show 
Jeopardy). Neural network models taking the form of a many-body potential (6.6.1) 
are referred to variously as higher-order Boltzmann machines [Sej86], higher-order 
neural networks [GM87] or simply spin-glass models [AGS85]. The general form of 
the quantity analogous to the potential in neural networks is given by 
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where T is the graph of the neural network, the Vj are weights and the s* take binary 
values [Sej86]. Models used in most applications are optimised in order to make 
decisions, and so are deterministic. However, by optimising the paramaters of (6.6.6) 
over a state space of inputs and outputs, s,, the optimisation phase of these models 
involves the artificial brain exploring a subset of the randomised phase space of the 
model. 
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